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Abstract 

These lecture notes in Lie Groups are designed for a 1-semester third year 
or graduate course in mathematics, physics, engineering, chemistry or biology. 
This landmark theory of the 20th Century mathematics and physics gives a 
rigorous foundation to modern dynamics, as well as field and gauge theories 
in physics, engineering and biomechanics. We give both physical and medical 
examples of Lie groups. The only necessary background for comprehensive 
reading of these notes are advanced calculus and linear algebra. 
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1 Preliminaries: Sets, Maps and Diagrams 
1.1 Sets 

Given a map (or, a function) f : A ^ B, the set A is called the domain of /, and 
denoted Dom/. The set B is called the codomain of /, and denoted Cod/. The 
codomain is not to be confused with the range of f{A), which is in general only a 
subset of B (see [SJE]). 

A map / : X — )• y is called injective, or 1-1, or an injection, iff for every y in the 
codomain Y there is at most one x in the domain X with f{x) = y. Put another way, 
given X and x' in X, if f{x) = f{x'), then it follows that x = x'. A map f : X ^ Y 
is called surjective, or onto, or a surjection, iff for every y in the codomain Cod/ 
there is at least one x in the domain X with /(x) = y. Put another way, the range 
f{X) is equal to the codomain Y . A map is bijective iff it is both injective and 
surjective. Injective functions are called monomorphisms, and surjective functions 
are called epimorphisms in the category of sets (see below). Bijective functions are 
called isomorphisms. 

A relation is any subset of a Cartesian product (see below). By definition, an 
equivalence relation a on a set X is a relation which is reflexive, symmetrical and 
transitive, i.e., relation that satisfies the following three conditions: 

1. Reflexivity: each element x G X is equivalent to itself, i.e., xax; 

2. Symmetry: for any two elements a,b G X, aab implies baa; and 

3. Transitivity: aab and bac implies aac. 

Similarly, a relation < defines a partial order on a set 5* if it has the following 
properties: 

1. Reflexivity: a < a for all a £ S; 

2. Antisymmetry: a <b and b < a implies a = b; and 

3. Transitivity: a < b and b < c implies a < c. 

A partially ordered set (or poset) is a set taken together with a partial order on 
it. Formally, a partially ordered set is defined as an ordered pair P = [X, <), where 
X is called the ground set of P and < is the partial order of P. 
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1.2 Maps 

Let / and g be maps with domains A and B. Then the maps f + g, f — g, fg, and 
f/g are defined as follows (see [8119]): 

[f + g){x) = f{x) + g{x) domain =Af^B, 

if ~ 9){^) = f{^) ~ 9{x) domain = An B, 

{f9){x) = f{x)g{x) domain =Ar\B, 

f\ fix) 

^ 'x) = —-4- domain = {x e Ar\B : g{x) / 0}. 



9 J 9[x) 

Given two maps / and g, the composite map / o <?, called the composition of / 
and g, is defined by 

if o 9){x) = f{9{x)). 

The (/o (^)— machine is composed of the (7— machine (first) and then the /—machine, 

x^[[g]]^g{x)^[[f]]^f{g{x)). 

For example, suppose that y = f{u) = y/u and u = g{x) = + 1. Since y is a 
function of u and u is a function of x, it follows that y is ultimately a function of x. 
We calculate this by substitution 



y 



f(u) = fog = f{g{x)) = f{x^ + 1) = Vx^ + 1. 



If / and g are both differentiable (or smooth, i.e., C°°) maps and h = fog is the 
composite map defined by h(x) = f{g{x)), then h is differentiable and h' is given by 
the product: 

h'{x) = ng{x))g'{x). 

In Leibniz notation, if y = f{u) and u = g{x) are both differentiable maps, then 

dy dy du 
dx du dx 

The reason for the name chain rule becomes clear if we add another link to the chain. 
Suppose that we have one more differentiable map x = h{t). Then, to calculate the 
derivative of y with respect to t, we use the chain rule twice, 

dy dy du dx 
dt du dx dt 



tan 



d{x,...) 

maps a region S onto a region R, we have the following substitution formulas: 



Given a 1-1 continuous (i.e., C^) map F with a nonzero Jacobi 



that 
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1. For a single integral. 



f dx 
f{x)dx= / f{x{u))—du; 
R Js du 



2. For a double integral, 



R 



f{x,y)dA= // f{x{u,v),y{u,v)) 
J Js 

3. For a triple integral, 
f{x,y,z)dV = 



d{x,y) 



d{u, v) 




f{x{u, V, w),y{u, V, w),z{u, v, w)) 



dudv; 



d{x,y,z) 



d{u, V, w) 



dudvdw: 



IR JJJs 
4. Generalization to n— tuple integrals is obvious. 

1.3 Commutative Diagrams 

Many properties of mathematical systems can be unified and simplified by a presen- 
tation with commutative diagrams of arrows. Each arrow f : X ^ Y represents a 
function (i.e., a map, transformation, operator); that is, a source (domain) set X, a 
target (codomain) set Y, and a rule x i— )• f{x) which assigns to each element x £ X 
an element f{x) gY. A typical diagram of sets and functions is (see [8l|9]): 



X 



f 



Y 



X 



f 



■fix) 




or 




This diagram is commutative iff /i = gof, where gofis the usual composite function 
g o f : X ^ Z, defined by x i— )• g{f{x)). 

Less formally, composing maps is like following directed paths from one object to 
another (e.g., from set to set). In general, a diagram is commutative iff any two paths 
along arrows that start at the same point and finish at the same point yield the same 
'homomorphism' via compositions along successive arrows. Commutativity of the 
whole diagram follows from commutativity of its triangular components. Study of 
commutative diagrams is popularly called 'diagram chasing', and provides a powerful 
tool for mathematical thought. 

Many properties of mathematical constructions may be represented by universal 
properties of diagrams. Consider the Cartesian product X xY oi two sets, consisting 
as usual of all ordered pairs (x, y) of elements x £ X and y € Y. The projections 
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{x,y) I—)- X, (x,y) I— 7- y of the product on its 'axes' X and Y are functions p : 
X xY ^ X, q : X xY ^ Y . Any function h : W ^ X x Y from a third set W 
is uniquely determined by its composites po h and q o h. Conversely, given W and 
two functions / and g as in the diagram below, there is a unique function h which 
makes the following diagram commute: 




This property describes the Cartesian product X x Y uniquely; the same diagram, 
read in the category of topological spaces or of groups, describes uniquely the Carte- 
sian product of spaces or of the direct product of groups. 



2 Groups 

A group is a pointed set {G, e) with a multiplication fi : G x G ^ G and an inverse 
V : G ^ G such that the following diagrams commute (see [lOl El E] ) : 



1. 



(e, 1) (l,e) 
G^^GxG^^^G 



I 

G 



(e is a two-sided identity) 



2. 



GxGx G ' Gx G 



1 X fl 

GxG 



G 



(associativity) 
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3. 




G 



(inverse). 

Here e : G — ^ G is the constant map e{g) = e for all g & G. (e, 1) means the map 
such that (e, l){g) = ie,g), etc. A group G is called commutative or Abelian group 
if in addition the following diagram commutes 

G X G -G X G 




G 



where T: GxG^GxG is the switch map T{gi,g2) = (52,51), for all (51,52) € 
G X G. 

A group G acis (on the left) on a set A if there is a function a : G x ^ — )• A 
such that the following diagrams commute: 



1. 




2. 



GxGxA . GxA 



fix 1 



GxA- 



a 



a 



A 



where (e, l)(x) = {e,x) for all x & A. The orbits of the action are the sets 
Gx = {gx : 5 G G} for all x e A. 
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Given two groups (G, *) and (H, •), a group homomorphism from (G, *) to {H, •) 
is a function h : G ^ H such that for all x and y in G it holds that 

h{x *y) = h{x) ■ h{y). 

Prom this property, one can deduce that h maps the identity element ec of G to 
the identity element en of H, and it also maps inverses to inverses in the sense that 
h{x~^) = h{x)~^. Hence one can say that h is compatible with the group structure. 

The kernel Ker /i of a group homomorphism h : G ^ H consists of all those 
elements of G which are sent by h to the identity element of H, i.e., 

Ker h = {x £ G : h{x) = en}- 

The image Imh of a group homomorphism h : G ^ H consists of all elements 
of G which are sent by h to i.e., 

Im/i = {h{x) : X G G}. 

The kernel is a normal subgroup of G and the image is a subgroup of H. The 
homomorphism h is injective (and called a group monomorphism) iff Ker/i = ec, 
i.e., iff the kernel of h consists of the identity element of G only. 

3 Manifolds 

A manifold is an abstract mathematical space, which locally (i.e., in a close-up 
view) resembles the spaces described by Euclidean geometry, but which globally 
(i.e., when viewed as a whole) may have a more complicated structure (see |11]). 
For example, the surface of Earth is a manifold; locally it seems to be flat, but viewed 
as a whole from the outer space (globally) it is actually round. A manifold can be 
constructed by 'gluing' separate Euclidean spaces together; for example, a world 
map can be made by gluing many maps of local regions together, and accounting 
for the resulting distortions. 

As main pure-mathematical references for manifolds we recommend popular 
graduate textbooks by two ex~Bourbaki members. Serge Lang [131 H] and Jean 
Dieudonne [14^ 115]. Besides, the reader might wish to consult some other 'classics', 
including [TTl [16l [iTl [18l [19l [201 E]- Finally, as first-order applications, we recom- 
mend three popular textbooks in mechanics, [2l[Tl|4], as well as our own geometrical 
monographs [H |9] . 

Another example of a manifold is a circle S^. A small piece of a circle appears to 
be like a slightly-bent part of a straight line segment, but overall the circle and the 
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segment are different ID manifoids. A circle can be formed by bending a straiglit 
line segment and gluing the ends together^ 

The surfaces of a spher^ and a torwjf] are examples of 2D manifolds. Manifolds 

^Locally, the circle looks like a line. It is ID, that is, only one coordinate is needed to say 
where a point is on the circle locally. Consider, for instance, the top part of the circle, where the 

coordinate is positive. Any point in this part can be described by the x— coordinate. So, there is a 
continuous bijection Xtop mapping which is 1-1 both ways) , which maps the top part of the circle 
to the open interval ( — 1, 1), by simply projecting onto the first coordinate: Xtopi^yV) ~ ^- Such 
a function is called a chart. Similarly, there are charts for the bottom, left , and right parts of the 
circle. Together, these parts cover the whole circle and the four charts form an atlas (see the next 
subsection) for the circle. The top and right charts overlap: their intersection lies in the quarter 
of the circle where both the x~ and the coordinates are positive. The two charts Xtop ^^i^ 
Xright map thls part bijectively to the interval (0, 1). Thus a function T from (0, 1) to itself can be 
constructed, which first inverts the top chart to reach the circle and then follows the right chart 
back to the interval: 

T{a) = Xright {Xtop{a)) = Xright («. - a^) = \/l-a^- 

Such a function is called a transition map. The top, bottom, left, and right charts show that the 
circle is a manifold, but they do not form the only possible atlas. Charts need not be geometric 
projections, and the number of charts is a matter of choice. T and the other transition functions 
are differentiable on the interval (0, 1). Therefore, with this atlas the circle is a differentiable, or 
smooth manifold. 

^The surface of the sphere can be treated in almost the same way as the circle . It can be 
viewed as a subset of R^, defined by: S = {{x, y, z) £ R^|a;^ + + = 1}. The sphere is 2D, so 
each chart will map part of the sphere to an open subset of M^. Consider the northern hemisphere, 
which is the part with positive z coordinate. The function x defined by xi^illyz) = {x,y), maps 
the northern hemisphere to the open unit disc by projecting it on the {x, y)— plane. A similar chart 
exists for the southern hemisphere. Together with two charts projecting on the (x, z)— plane and two 
charts projecting on the (y, z)— plane, an atlas of six charts is obtained which covers the entire sphere. 
This can be easily generalized to an nD sphere S" = {(a;i, 2:2, a;„) G R"|a;i + X2 + ... + x^ — 1}. 

An n— sphere 5" can be also constructed by gluing together two copies of M". The transition 
map between them is defined as R" \ {0} — )■ R" \ {0} : x 1— >■ a:;/||a;|p. This function is its own inverse, 
so it can be used in both directions. As the transition map is a {C°°)~ smooth function, this atlas 
defines a smooth manifold. 

torus (pi. tori), denoted by T^, is a doughnut-shaped surface of revolution generated by 
revolving a circle about an axis coplanar with the circle. The sphere is a special case of the 
torus obtained when the axis of rotation is a diameter of the circle. If the axis of rotation does not 
intersect the circle, the torus has a hole in the middle and resembles a ring doughnut, a hula hoop 
or an infiated tire. The other case, when the axis of rotation is a chord of the circle, produces a 
sort of squashed sphere resembling a round cushion. 

A torus can be defined parametrically by: 

x{u,v) — {R + r cosv) cosu, y{u,v) = {R + rcosv)sinu, z{u,v) = rsinv, 

where u,v £ [0, 2tt], R is the distance from the center of the tube to the center of the torus, and r 
is the radius of the tube. According to a broader definition, the generator of a torus need not be a 
circle but could also be an ellipse or any other conic section. 

Topologically, a torus is a closed surface defined as product of two circles: — x . The 
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are important objects in mathematics, physics and control theory, because they 
allow more complicated structures to be expressed and understood in terms of the 
well-understood properties of simpler Euclidean spaces (see [9]). 

The Cartesian product of manifolds is also a manifold (note that not every man- 
ifold can be written as a product). The dimension of the product manifold is the 
sum of the dimensions of its factors. Its topology is the product topology, and a 
Cartesian product of charts is a chart for the product manifold. Thus, an atlas 
for the product manifold can be constructed using atlases for its factors. If these 
atlases define a differential structure on the factors, the corresponding atlas defines 
a differential structure on the product manifold. The same is true for any other 
structure defined on the factors. If one of the factors has a boundary, the product 
manifold also has a boundary. Cartesian products may be used to construct tori 
and cylinders, for example, as x and x [0, 1], respectively. 

Manifolds need not be connected (all in 'one piece'): a pair of separate circles 
is also a topological manifold{see below). Manifolds need not be closed: a line 
segment without its ends is a manifold. Manifolds need not be finite: a parabola is 
a topological manifold. 

Manifold^ can be viewed using either extrinsic or intrinsic view. In the extrinsic 
view, usually used in geometry and topology of surfaces, an nD manifold M is seen 
as embedded in an (n + 1)D Euclidean space M""*"^. Such a manifold is called a 
'codimension 1 space'. With this view it is easy to use intuition from Euclidean 
spaces to define additional structure. For example, in a Euclidean space it is always 
clear whether a vector at some point is tangential or normal to some surface through 
that point. On the other hand, the intrinsic view of an nD manifold M is an 
abstract way of considering M as a topological space by itself, without any need 

surface described above, given the relative topology from R'^, is homeomorphic to a topological torus 
as long as it does not intersect its own axis. 

One can easily generalize the torus to arbitrary dimensions. An n— torus T" is defined as a 
product of n circles: T" — x x ■ ■ ■ x . Equivalently, the n— torus is obtained from the 
n— cube (the R"— generalization of the ordinary cube in R^) by gluing the opposite faces together. 

An n— torus T" is an example of an nD compact manifold. It is also an important example of a 
Lie group (see below). 

4 

Additional structures are often defined on manifolds. Examples of manifolds with additional 
structure include: 

• dijferentiable (or, smooth manifolds, on which one can do calculus; 

• Riemannian manifolds, on which distances and angles can be defined; they serve as the 
configuration spaces in mechanics; 

• symplectic manifolds, which serve as the phase spaces in mechanics and physics; 

• 4D pseudo-Riemannian manifolds which model space-time in general relativity. 
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for surrounding (n + 1)D Euclidean space. This view is more flexible and thus it 
is usually used in high-dimensional mechanics and physics (where manifolds used 
represent configuration and phase spaces of dynamical systems) , can make it harder 
to imagine what a tangent vector might be. 

3.1 Definition of a Manifold 

Consider a set M (see Figure [1]) which is a candidate for a manifold. Any point 
x S M has its Euclidean chart, given by a 1-1 and onto map ip^ : M ^ M", with its 
Euclidean image Vi = ip^{U,j). More precisely, a chart (p^ is defined by (see [SlIS]) 

93i : M D [/, 9 X V?,(x) G C M", 

where Ui <Z M and Vi C M" are open sets. 




Figure 1: Geometric picture of the manifold concept. 

Clearly, any point x G M can have several different charts (see Figure [T]). Con- 
sider a case of two charts, v'i) ■ — ?• M", having in their images two open sets, 
Vij = iPi{Ui n Uj) and Vji = ^j{Ui n Uj). Then we have transition functions ip^j 
between them, 

'^ij = "Pj ° ■ Vij yji, locally given by (fij{x) = (pj{(pY^{x)). 

If transition functions ipj_j exist, then we say that two charts, ip^ and ipj are compat- 
ible. Transition functions represent a general (nonlinear) transformations of coordi- 
nates, which are the core of classical tensor calculus. 
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A set of compatible charts ip^ : M ^ M", such that each point x £ M has its 
Euclidean image in at least one chart, is called an atlas. Two atlases are equivalent 
iff all their charts are compatible (i.e., transition functions exist between them), so 
their union is also an atlas. A manifold structure is a class of equivalent atlases. 

Finally, as charts ip^ : M ^ were supposed to be 1-1 and onto maps, they 
can be either homeomorphisms, in which case we have a topological (C^) manifold, 
or diffeomorphisms, in which case we have a smooth (C^) manifold. 

3.2 Formal Definition of a Smooth Manifold 

Given a chart {U,(p), we call the set U a coordinate domain, or a coordinate neigh- 
borhood of each of its points. If in addition (p{U) is an open ball in M", then U 
is called a coordinate ball. The map if is called a (local) coordinate map, and the 
component functions {x^ , x^) of (f, defined by ip{m) = (x^(m), a:"(m)), are 
called local coordinates on f7 [HI [9] . 

Two charts {Ui,ipi) and {U2,f2) such that C/i H C/2 7^ are called compatible 
if ^i{Ui n U2) and 932(^2 H Ui) are open subsets of M". A family {Ua,^a)aeA of 
compatible charts on M such that the Ua form a covering of M is called an atlas. 
The maps = ipp o ip~^ : Lp^{Uap) ipp{Uap) are called the transition maps, 
for the atlas {Ua,^a)a&A-, where Uap = Ua ri Up, so that we have a commutative 
triangle: 

Uap C M 

'PaiUal^) fpiUap) 

An atlas {Ua,'Pa)a£A for a manifold Af is said to be a C^— atlas, if all transi- 
tion maps ipai3 '■ (faiUafs) — ^ ^i^iUa/s) are of class C*^. Two atlases are called 
C^ — equivalent, if their union is again a C*^— atlas for M. An equivalence class of 
C'^— atlases is called a C^— structure on M. In other words, a smooth structure on 
M is a maximal smooth atlas on M, i.e., such an atlas that is not contained in any 
strictly larger smooth atlas. By a C''— manifold M, we mean a topological manifold 
together with a C'^— structure and a chart on M will be a chart belonging to some 
atlas of the C*^— structure. Smooth manifold means C°°— manifold, and the word 
^ smooth^ is used synonymously with C°°. 

Sometimes the terms 'local coordinate system' or 'parametrization' are used 
instead of charts. That M is not defined with any particular atlas, but with an 
equivalence class of atlases, is a mathematical formulation of the general covariance 
principle. Every suitable coordinate system is equally good. A Euclidean chart 
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may well suffice for an open subset of M", but this coordinate system is not to be 
preferred to the others, which may require many charts (as with polar coordinates), 
but are more convenient in other respects. 

For example, the atlas of an n— sphere 5" has two charts. If N = (1,0, ...,0) 
and S" = (— 1, 0, 0) are the north and south poles of S'"' respectively, then the two 
charts are given by the stereographic projections from N and S: 



V^i : 5"\{Af}^M",(pi(x\...,2;"+i) = (xV(l-2;^),---,2;"+V(l-xi)), and 
: S^\{S}^R^,ip,ix\...,x^-^') = ix^/il + x'),...,x^+'/il+x')), 



while the overlap map ^2 ° '■ R"'\{0} — )• M"\{0} is given by the diffeomorphism 
if 2 ° fi^){z) = -2/1 kiP, for z in M"\{0}, from M"\{0} to itself. 

Various additional structures can be imposed on R", and the corresponding man- 
ifold M will inherit them through its covering by charts. For example, if a covering 
by charts takes their values in a Banach space E, then E is called the model space 
and M is referred to as a C''— Banach manifold modelled on E. Similarly, if a cov- 
ering by charts takes their values in a Hilbert space Ti, then T-L is called the model 
space and M is referred to as a C''— Hilbert manifold modelled on T-L. If not other- 
wise specified, we will consider M to be an Euclidean manifold, with its covering by 
charts taking their values in R". 

For a Hausdorff C'^— manifold the following properties are equivalent: (i) it is 
paracompact; (ii) it is metrizable; (iii) it admits a Riemannian metric jl (iv) each 
connected component is separable. 

3.3 Smooth Maps Between Smooth Manifolds 

A map if : M ^ N between two manifolds M and A^, with M 9 m 1— )• ip{m) E N, is 
called a smooth map, or C'^— map, if we have the following charting [8l|9]: 

^Recall the corresponding properties of a Euclidean metric d. For any three points x,y,z £ R", 
the following axioms are valid: 



Ml 
Ma 



d{x,y) > 0, for x / y; 
d{x,y) = d{y,x); 



and d{x, y) = 0, for x — y; 
M3 : d{x,y) < d{x,z) + d{z,y). 
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\ M 


•P 


N f 










(j){m) 



V 



1p O If O (p 



-1 



0, 




This means that for each m €z M and each chart {V, ip) on N with ip (m) G V there 
is a chart {U, (j)) on M with m £ U,ip {U) C T/, and $ 
the fohowing diagram commutes: 



'ip o (f o (f) ^ is C"^, that is, 



V CN 



Let M and be smooth manifolds and let (/9 : M — t- be a smooth map. The 
map (f is called a covering, or equivalently, M is said to cover N, if (/9 is surjective 
and each point n £ N admits an open neighborhood V such that ip~^{V) is a union 
of disjoint open sets, each diffeomorphic via (p to V. 

A C^— map if : AI ^ N is called a — diffeomorphism if is a bijection, 
: N ^ M exists and is also C^. Two manifolds are called diffeomorphic if there 
exists a diffeomorphism between them. All smooth manifolds and smooth maps 
between them form the category 7W. 



3.4 Tangent Bundle and Lagrangian Dynamics 

The tangent bundle of a smooth n— manifold is the place where tangent vectors live, 
and is itself a smooth 2n— manifold. Vector-fields are cross-sections of the tangent 
bundle. The Lagrangian is a natural energy function on the tangent bundle (see 

mm)- 
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In mechanics, to each nD configuration manifold M there is associated its 2nD 
velocity phase-space manifold, denoted by TM and cahed the tangent bundle of 
M (see Figure [2]). The original smooth manifold M is called the base of TM. 
There is an onto map tt : TM — )• M, called the projection. Above each point 
X £ M there is a tangent space T^M = 'ir~^{x) to M at x, which is called a fibre. 
The fibre T^M C TM is the subset of TM, such that the total tangent bundle, 
TM = I I TxM, is a disjoint union of tangent spaces TxM to M for all points 

mGA/ 

X G M. From dynamical perspective, the most important quantity in the tangent 
bundle concept is the smooth map u : M — )• TM, which is an inverse to the projection 
TT, i.e, vr o ti = Mm, vr(ti(x)) = x. It is called the velocity vector-field. Its graph 
{x,v{x)) represents the cross-section of the tangent bundle TM. This explains 
the dynamical term velocity phase-space, given to the tangent bundle TM of the 
manifold M. 




M 



X 



Figure 2: A sketch of a tangent bundle TM of a smooth manifold M. 

If [a, b] is a closed interval, a C"— map 7 : [a, 6] — )■ M is said to be differentiable 
at the endpoint a if there is a chart (U, (p) at 7(a) such that the following limit exists 
and is finite: 



^((/> o 7) (a) = {4)0 7)' (a) = lim ■ 
at t^a 



- (0o7)(q) 
t — a 



(1) 



Generalizing ([T]), we get the notion of the curve on a manifold. For a smooth 
manifold M and a point m G Af a curve at m is a C''— map 7 : / — )• M from an 
interval / C M into M with G / and 7(0) = m. 

Two curves 7^ and 72 passing though a point m G U are tangent at m with 
respect to the chart (C/, if (</'O7i)'(0) = (</'O72)'(0)- Thus, two curves are tangent 
if they have identical tangent vectors (same direction and speed) in a local chart on 
a manifold. 

For a smooth manifold M and a point m G M, the tangent space TmM to M at 
m is the set of equivalence classes of curves at m: 

TmM = {[y]m : 7 is a curve at a point m G M}. 
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A C'^— map if : M 3 € N between two manifolds M and N induces 

a linear map Tjnip : T^M — >■ T^^jj^-jN for each point m G M, called a tangent map, 
if we have: 

T^{M) TM 



TTM 



T{N) T^^ra){N) 



©■ 



m 



(5 



71" AT 




i.e., the following diagram commutes: 



TraM. 



Tm'P 



M 3 m- 



TTN 



(p{m) e N 



with the natural projection ttm '■ TM — > M, given by TTMiTmM) = m, that takes a 
tangent vector v to the point m G M at which the vector v is attached i.e., v G Tj^M. 

For an nD smooth manifold M, its nD tangent bundle TM is the disjoint union 
of all its tangent spaces T^M at all points m G M, TM = j | T^M. 

meM 

To define the smooth structure on TM, we need to specify how to construct 
local coordinates on TM. To do this, let x"'{m)) be local coordinates of a 

point m on M and let {v^{m), ...,f"(m)) be components of a tangent vector in this 
coordinate system. Then the 2n numbers {x^{m), ...,x'^{m), v^{m), ...,v'^{m)) give 
a local coordinate system on TM. 

TM = |_J TmM defines a family of vector spaces parameterized by M. The 

mGM 

inverse image TT~^{m) of a point m € M under the natural projection ttm is the 
tangent space T^M. This space is called the fibre of the tangent bundle over the 
point m G M. 
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A C'^-map (/? : M 
tangent map Tip : TM 
diagram commutes: 



between two manifolds M and N induces a linear 
TN between their tangent bundles, i.e., the following 



TM 



TN 




TTN 



All tangent bundles and their tangent maps form the category TB. The category 
TB is the natural framework for Lagrangian dynamics. 

Now, we can formulate the global version of the chain rule. If (/? : M — t- A and 
ip : N ^ P are two smooth maps, then we have T{tp o ip) = Tip o Tip. In other 
words, we have a functor T : Ai ^ TB from the category M of smooth manifolds 
to the category TB of their tangent bundles: 

TM 



TN 





3.5 Cotangent Bundle and Hamiltonian Dynamics 

The cotangent bundle of a smooth ?7,— manifold is the place is where 1-forms live, and 
is itself a smooth 2n— manifold. Covector-fields (l~forms) are cross-sections of the 
cotangent bundle. The Hamiltonian is a natural energy function on the cotangent 
bundle (see [8l[9]). 

A dual notion to the tangent space T^M to a smooth manifold M at a point m 
is its cotangent space T^M at the same point m. Similarly to the tangent bundle, 
for a smooth manifold M of dimension n, its cotangent bundle T*M is the disjoint 
union of all its cotangent spaces T^M at all points m £ M, i.e., T*M = | | T^M. 

Therefore, the cotangent bundle of an n— manifold M is the vector bundle T*M = 
(TM)*, the (real) dual of the tangent bundle TM. 

If M is an n— manifold, then T*M is a 2n— manifold. To define the smooth 
structure on T*M, we need to specify how to construct local coordinates on T*M. 
To do this, let (x^(m), a;"(m)) be local coordinates of a point m on M and let 
(pi(m), ...,pn{m)) be components of a covector in this coordinate system. Then the 
2n numbers {x^{m), x'^{m), pi(m), ...,p„(m)) give a local coordinate system on 
T*M. This is the basic idea one uses to prove that indeed T*M is a 2n— manifold. 
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T*M = I I T^M defines a family of vector spaces parameterized by M, with 

meM 

tlie conatural projection, vr^^ : T*M — )• M, given by 7r^(T^M) = m, that takes a 
covector p to the point m £ M at which the covector p is attached i.e., p S T^M. 
The inverse image -K^^im) of a point m S M under the conatural projection vr^^ is 
the cotangent space T^M. This space is called the fibre of the cotangent bundle 
over the point m & M. 

In a similar way, a C'^— map ip : M ^ N between two manifolds M and induces 
a linear cotangent map T*ip : T*M — )■ r*A^ between their cotangent bundles, i.e., 
the following diagram commutes: 



T*M 



T*Lp 



T*N 



M 



N 



M 



N 



All cotangent bundles and their cotangent maps form the category T*B. The 
category T*B is the natural stage for Hamiltonian dynamics. 

Now, we can formulate the dual version of the global chain rule, li ip : M ^ N 
and ip : N ^ P are two smooth maps, then we have T*(^ o ip) = T*'ip o T*ip. In 
other words, we have a cofunctor T* : 7W =^ T*B from the category Ai of smooth 
manifolds to the category T*B of their cotangent bundles: 

T*M 





4 Lie Groups 

In this section we present the basics of classical theory of Lie groups and their Lie 
algebras, as developed mainly by Sophus Lie, Elie Cartan, Felix Klein, Wilhelm 
Killing and Hermann Weyl. For more comprehensive treatment see e.g., |2 H 122 ^ [23 | 
[22 [25]. 

In the middle of the 19th Century S. Lie made a far reaching discovery that 
techniques designed to solve particular unrelated types of ODEs, such as separable, 
homogeneous and exact equations, were in fact all special cases of a general form of 
integration procedure based on the invariance of the differential equation under a 
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continuous group of symmetries. Roughly speaking a symmetry group of a system 
of differential equations is a group that transforms solutions of the system to other 
solutions. Once the symmetry group has been identified a number of techniques 
to solve and classify these differential equations becomes possible. In the classical 
framework of Lie, these groups were local groups and arose locally as groups of 
transformations on some Euclidean space. The passage from the local Lie group to 
the present day definition using manifolds was accomplished by E. Cartan at the 
end of the 19th Century, whose work is a striking synthesis of Lie theory, classical 
geometry, differential geometry and topology. 

These continuous groups, which originally appeared as symmetry groups of dif- 
ferential equations, have over the years had a profound impact on diverse areas such 
as algebraic topology, differential geometry, numerical analysis, control theory, clas- 
sical mechanics, quantum mechanics etc. They are now universally known as Lie 
groups. 

A Lie group is smooth manifold which also carries a group structure whose 
product and inversion operations are smooth as maps of manifolds. These objects 
arise naturally in describing physical symmetriesjf] 

^Here are a few examples of Lie groups and their relations to other areas of mathematics and 
physics: 

1. Euclidean space R" is an Abelian Lie group (with ordinary vector addition as the group 
operation). 

2. The group GL„(R) of invertible matrices (under matrix multiplication) is a Lie group of 
dimension . It has a subgroup SL„(R) of matrices of determinant 1 which is also a Lie 
group. 

3. The group On(R) generated by all rotations and reflections of an nD vector space is a Lie 
group called the orthogonal group. It has a subgroup of elements of determinant 1, called 
the special orthogonal group SO{n), which is the group of rotations in R". 

4. Spin groups are double covers of the special orthogonal groups (used e.g., for studying 
fermions in quantum field theory). 

5. The group Sp2n(R) of all matrices preserving a symplectic form is a Lie group called the 
symplectic group. 

6. The Lorentz group and the Poincare group of isometrics of space-time are Lie groups of 
dimensions 6 and 10 that are used in special relativity. 

7. The Heisenberg group is a Lie group of dimension 3, used in quantum mechanics. 

8. The unitary group U{n) is a compact group of dimension consisting of unitary matrices. 
It has a subgroup of elements of determinant 1, called the special unitary group SU{n). 

9. The group U{1) x SU{2) x 5(7(3) is a Lie group of dimension 1 + 3 + 8 = 12 that is the 
gauge group of the Standard Model of elementary particles, whose dimension corresponds to: 
1 photon + 3 vector bosons + 8 gluons. 
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A Lie group is a group whose elements can be continuously parametrized by 
real numbers, such as the rotation group 50(3), which can be parametrized by the 
Euler angles. More formally, a Lie group is an analytic real or complex manifold 
that is also a group, such that the group operations multiplication and inversion 
are analytic maps. Lie groups are important in mathematical analysis, physics and 
geometry because they serve to describe the symmetry of analytical structures. They 
were introduced by Sophus Lie in 1870 in order to study symmetries of differential 
equations. 

While the Euclidean space is a real Lie group (with ordinary vector addition 
as the group operation), more typical examples are given by matrix Lie groups, i.e., 
groups of invertible matrices (under matrix multiplication). For instance, the group 
S0{3) of all rotations in is a matrix Lie group. 

One classifies Lie groups regarding their algebraic propertied (simple, semisim- 
ple, solvable, nilpotent, Abelian), their connectedness (connected or simply con- 
nected) and their compactnesJ^ 

To every Lie group, we can associate a Lie algebra which completely captures 
the local structure of the group (at least if the Lie group is connected) H 

^If G and H are Lie groups (both real or both complex), then a Lie-group-homomorphism 
/ : G — > -ff is a group homomorphism which is also an analytic map (one can show that it is 
equivalent to require only that f be continuous). The composition of two such homomorphisms 
is again a homomorphism, and the class of all (real or complex) Lie groups, together with these 
morphisms, forms a category. The two Lie groups are called isomorphic iff there exists a bijective 
homomorphism between them whose inverse is also a homomorphism. Isomorphic Lie groups do 
not need to be distinguished for all practical purposes; they only differ in the notation of their 
elements. 

*An n— torus T" = x x ■ ■ ■ x (as defined above) is an example of a compact Abelian 
Lie group. This follows from the fact that the unit circle is a compact Abelian Lie group (when 
identified with the unit complex numbers with multiplication) . Group multiplication on T" is then 
defined by coordinate-wise multiplication. 

Toroidal groups play an important part in the theory of compact Lie groups. This is due in part 
to the fact that in any compact Lie group one can always find a maximal torus; that is, a closed 
subgroup which is a torus of the largest possible dimension. 

^Conventionally, one can regard any field X of tangent vectors on a Lie group as a partial 
differential operator, denoting by Xf the Lie derivative (the directional derivative) of the scalar 
field / in the direction of X. Then a vector-field on a Lie group G is said to be left-invariant if it 
commutes with left translation, which means the following. Define Lg[f]{x) — f{gx) for any analytic 
function / : G — >■ R and all g,x £ G. Then the vector-field X is left-invariant iff XLg = LgX for 
all g £ G. Similarly, instead of R, we can use C. The set of all vector-fields on an analytic manifold 
is a Lie algebra over R (or C). 

On a Lie group G, the left-invariant vector-fields form a subalgebra, the Lie algebra g associated 
with G. This Lie algebra is finite-dimensional (it has the same dimension as the manifold G) which 
makes it susceptible to classification attempts. By classifying g, one can also get a handle on the 
group G. The representation theory of simple Lie groups is the best and most important example. 

Every element v of the tangent space Te at the identity element e of G determines a unique 
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4.1 Definition of a Lie Group 



A Lie group is a smooth (Banach) manifold M that has at the same time a group 
G— structure consistent with its manifold M— structure in the sense that group mul- 
tiplication fi : G X G ^ G, {g,h) gh and the group inversion v : G ^ G, g ^ 
g~^ are C'^— maps. A point e G G is called the group identity element (see e.g., 

For example, any nD Banach vector space V is an Abelian Lie group with group 
operations fi : V xV ^ V, fi{x, y) = x + y, and u : V ^ V, u{x) = —x. The identity 
is just the zero vector. We call such a Lie group a vector group. 

Let G and H be two Lie groups. A map G — )• is said to be a morphism of 
Lie groups (or their smooth homomorphism) if it is their homomorphism as abstract 
groups and their smooth map as manifolds. 

Similarly, a group G which is at the same time a topological space is said to be 
a topological group if both maps (jU, z^) are continuous, i.e., C'^— maps for it. The 
homomorphism G ^ H of topological groups is said to be continuous if it is a 
continuous map. 

A topological group (as well as a smooth manifold) is not necessarily Hausdorff. 

left -invariant vector-field whose value at the element g of G is denoted by gv; the vector space 
underlying the Lie algebra g may therefore be identified with Te- 

Every vector-field v in the Lie algebra g determines a function c : R — > G whose derivative 
everywhere is given by the corresponding left-invariant vector-field: c'(t) — TL^-fffV and which 
has the property: c(s + t) — c{s)c{t), (for all s and t) (the operation on the r.h.s. is the group 
multiplication in G). The formal similarity of this formula with the one valid for the elementary 
exponential function justifies the definition: exp(tj) — c(l). This is called the exponential map, 
and it maps the Lie algebra g into the Lie group G. It provides a diffeomorphism between a 
neighborhood of in g and a neighborhood of e in G. This exponential map is a generalization of 
the exponential function for real numbers (since R is the Lie algebra of the Lie group of positive 
real numbers with multiplication), for complex numbers (since C is the Lie algebra of the Lie 
group of non-zero complex numbers with multiplication) and for matrices (since M(n, R) with the 
regular commutator is the Lie algebra of the Lie group GL(n, R) of all invertible matrices). As the 
exponential map is surjective on some neighborhood A'^ of e, it is common to call elements of the 
Lie algebra infinitesimal generators of the group G. 

The exponential map and the Lie algebra determine the local group structure of every connected 
Lie group, because of the Baker-Camphell-Hausdorff formula: there exists a neighborhood U of 
the zero element of the Lie algebra g, such that for u,v £ U we have 

exp(u)exp(ii) — exp(u + v + l/2[u, v] + l/12[[u, v],v] — 1/12[[m, v],u] — ...), 

where the omitted terms are known and involve Lie brackets of four or more elements. In case u 
and V commute, this formula reduces to the familiar exponential law: 

exp(«)exp(t;) = exp(it + v). 

Every homomorphism f : G H oi Lie groups induces a homomorphism between the corre- 
sponding Lie algebras g and f). The association G =^ g is called the Lie Functor. 
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A topological group G is HausdorfF iff its identity is closed. As a corollary we have 
that every Lie group is a HausdorfF topological group. 
For every 5 in a Lie group G, the two maps, 



^9 



G 
G 



G, 
G, 



h !->■ gh, 

9 1-^ gh, 



are called left and right translation maps. Since Lg o 



Lgh, and Rg o R,^ 



it follows that (Lg) 
phisms. Moreover Lg 



Lg-i and {Rg) = Rg 



so both Lg and Rg 



Rgh, 

are diffeomor- 



o Ru 



Rh o Lg, i.e. 



left and right translation commute. 



A vector-field X on G is called left-invariant vector-field if for every g ^ G, 
L*X = X, that is, if {ThLg)X{h) = X{gh) for all h e G, i.e., the following diagram 
commutes: 




A Riemannian metric on a Lie group G is called left-invariant if it is preserved by 
all left translations Lg, i.e., if the derivative of left translation carries every vector to 
a vector of the same length. Similarly, a vector field X on G is called left-invariant 
if (for every g e G) L*X = X. 



4.2 Lie Algebra 

An algebra A is a vector space with a product. The product must have the property 
that 

a{uv) = {au)v = u{av), 

for every a G M and u,v € A. A map (p : A ^ A' between algebras is called an 
algebra homomorphism, if (j){u ■ v) = (j){u) ■ <^(v). A vector subspace 3 of an algebra 
A is called a left ideal (resp. right ideal) if it is closed under algebra multiplication 
and if u e A and i e 3 implies that ui e (resp. iu e 3). A subspace J is said 
to be a two-sided ideal if it is both a left and right ideal. An ideal may not be an 
algebra itself, but the quotient of an algebra by a two-sided ideal inherits an algebra 
structure from A. 

A Lie algebra is an algebra A where the multiplication, i.e., the Lie bracket 
{u,v) I-)- [u,v], has the following properties: 

LA 1. [u, u] = for every u £ A, and 

LA 2. [u, [v, w]] + [w, [u, v]] + [v, w, u]\ = for all u,v,w G A. 
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The condition LA 2 is usually called Jacobi identity. A subspace E C A of a Lie 
algebra is called a Lie subalgebra if [u, v] & E for every u,v £ E. A map (p : A ^ A' 
between Lie algebras is called a Lie algebra homomorphism if (j){[u, v]) = [(j){u), (j){v)] 
for each u,v G A. 

All Lie algebras (over a given field IK) and all smooth homomorphisms between 
them form the category CAC, which is itself a complete subcategory of the category 
AC of all algebras and their homomorphisms. 

Let Xl[G) denote the set of left-invariant vector-fields on G; it is a Lie subalge- 
bra of X{G), the set of ah vector-fields on G, since L*g[X, Y] = [L*gX, L*gY] = [X, Y], 
so the Lie bracket [X,Y] G Xl{G). 

Let e be the identity element of G. Then for each ^ on the tangent space T^G we 
define a vector-field X^ on G by X^{g) = TeLg{^). Xl{G) and TgG are isomorphic 
as vector spaces. Define the Lie bracket on T^G by [^,77] = [X^,X^](e) for all 
^, ?? G TgG. This makes T^G into a Lie algebra. Also, by construction, we have 
[X^,Xri] = X[^^j^y, this defines a bracket in TeG via left extension. The vector space 
TeG with the above algebra structure is called the Lie algebra of the Lie group G 
and is denoted q. 

For example, let ^ be a nD vector space. Then TgV ~ V and the left-invariant 
vector-field defined by ^ G Tf.V is the constant vector-field X^{r]) = ^, for all r] G V. 
The Lie algebra of y is itself. 

Since any two elements of an Abelian Lie group G commute, it follows that all 
adjoint operators Adg, g £ G, equal the identity. Therefore, the Lie algebra g is 
Abelian; that is, [C,r]] = for all C,r] £ g. 

For example, G = 50(3) is the group of rotations of 3D Euclidean space, i.e. 
the configuration space of a rigid body fixed at a point. A motion of the body is 
then described by a curve g = g{t) in the group S0{3). Its Lie algebra g = so(3) is 
the 3D vector space of angular velocities of all possible rotations. The commutator 
in this algebra is the usual vector (cross) product (see, e.g. [ll[3l[9]). 

A rotation velocity g of the rigid body (fixed at a point) is a tangent vector to 
the Lie group G = SO (3) at the point g £ G. To get the angular velocity, we must 
carry this vector to the tangent space TGe of the group at the identity, i.e. to its Lie 
algebra g = so (3). This can be done in two ways: by left and right translation, Lg 
and Rg. As a result, we obtain two different vector fields in the Lie algebra so(3) : 

ujc = Lg-i^g £ so(3) and = Rg-ug S so(3), 

which are called the 'angular velocity in the body' and the 'angular velocity in space,' 
respectively. 

The dual space g* to the Lie algebra g = so (3) is the space of angular momenta 
TT. The kinetic energy T of a body is determined by the vector field of angular 
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velocity in the body and does not depend on the position of the body in space. 
Therefore, kinetic energy gives a left-invariant Riemannian metric on the rotation 
group G = SO (3). 

4.3 One-Parameter Subgroup 

Let be a left-invariant vector-field on G corresponding to in g. Then there is 
a unique integral curve 7^ : R — )• G of starting at e, i.e., (see, e.g. [El [9]) 

7g(t)=X^(75(t)), 7c(0) = e 

7^(t) is a smooth one-parameter subgroup of G, i.e., 7^(t + s) = 7^(t) • 7g(s), since, 
as functions of t both sides equal 7^(s) at t = and both satisfy differential equation 
j{t) = (7g(t)) by left invariance of X^, so they are equal. Left invariance can be 
also used to show that 7^(t) is defined for all t G M. Moreover, if </> : M — )• G is a 
one-parameter subgroup of G, i.e., a smooth homomorphism of the additive group 
M into G, then (/> = 7^ with ^ = 0(0), since taking derivative at s = in the relation 

(j){t + s) = ^{t) ■ 0(s) gives </.(t) = {(Pit)) , 

so = 7^ since both equal e at t = 0. Therefore, all one-parameter subgroups of G 
are of the form 7g(i) for some ^ G 0. 

4.4 Exponential Map 

The map exp : g — )• G, given by (see, e.g. [3l[8l|9]): 

exp(^) = 7^(1), exp(0) = e 

is called the exponential map of the Lie algebra g of G into G. exp is a G^ — 
map, similar to the projection vr of tangent and cotangent bundles; exp is locally a 
diffeomorphism from a neighborhood of zero in g onto a neighborhood of e in G; if 
/ : G — )• -ff is a smooth homomorphism of Lie groups, then 

/ o expc = ex.pH oTef . 

Also, in this case 

exp(sC) = 7^(s). 

Indeed, for fixed s G M, the curve t 1— )• 7^(ts), which at t = passes through e, 
satisfies the differential equation 

^75 (ts) = sX^ (7g(ts)) = Xs^ (75 (ts)) . 
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Since 7<jg(t) satisfies the same differential equation and passes through e at t = 0, it 
follows that 7s^(t) = 7g(st)- Putting t = 1 induces exp(s^) = 7^(s). 

Hence exp maps the line in g onto the one-parameter subgroup 7^(s) of G, 
which is tangent to ^ at e. It follows from left invariance that the flow of X 
satisfies F^{g) = gex.p{s^). 

Globally, the exponential map exp is a natural operation, i.e., for any morphism 
: G ^ H oi Lie groups G and H and a Lie functor F, the following diagram 
commutes: 

F{G) F{H) 



exp 



exp 



G jp H 

Let Gi and G2 be Lie groups with Lie algebras gi and Q2- Then Gi x G2 is a 
Lie group with Lie algebra gi x g2, and the exponential map is given by: 

exp : gi X g2 Gi X G2, (^1,^2) ^ (expi(^i), exp2(C2)) • 

For example, in case of a nD vector space, or infinite-dimensional Banach space, 
the exponential map is the identity. 

The unit circle in the complex plane S"^ = {z G C : |z| = 1} is an Abelian Lie 
group under multiplication. The tangent space Tf>S^ is the imaginary axis, and we 
identify M with TeS^ hy t ^ 2TTit. With this identification, the exponential map 
exp : R — )• S*-*^ is given by exp(t) = e^''**. 

The nD torus T" = S^x---xS^ (n times) is an Abelian Lie group. The exponen- 
tial map exp : — t- is given by 

exp(ti,...,t„) = (e2-*\...,e2-*"). 

Since = R/Z, it fohows that T" = R"/Z", the projection R" T" being given 
by the exp map. 

4.5 Adjoint Representation 

For every g £ G, the map (see, e.g. [H El El [9] ) : 

Adg = Te [Rg-i oLg) : g ^ g 

is called the adjoint map, or adjoint operator associated with g. 
For each ^ G g and g £ G we have 



exp (AdgS,) = (/(expOfi- 



-1 
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The relation between the adjoint map and the Lie bracket is the following: For 
all S,,r] £ Q we have 

t=o 



dt 



Left and right translations induce operators on the cotangent space T*Gg dual 
to Lg^ and denoted by (for every h £ G): 

g ■ -f (-Jgh — ^ J <^h, iig ■ J (-^/ig — ^ J L^/i- 

The transpose operators Ad*g : g — > g satisfy the relations Ad*f^ = Ad^^Ad* (for every 
g,h £ G) and constitute the co-adjoint representation of the Lie group G. The 
co-adjoint representation plays an important role in all questions related to (left) 
invariant metrics on the Lie group. According to A. Kirillov, the orbit of any vector 
field X in a Lie algebra g in a co-adjoint representation Ad* is itself a symplectic 
manifold and therefore a phase space for a Hamiltonian mechanical system. 

A Lie subgroup of G is a subgroup H of G which is also a submanifold of 
G. Then f) is a Lie subalgebra of g and moreover i) = £ q\ exp(t^) £ H, for all 
t £ R}. 

One can characterize Lebesgue measure up to a multiplicative constant on R" by 
its invariance under translations. Similarly, on a locally compact group there is a 
unique (up to a nonzero multiplicative constant) left-invariant measure, called Haar 
measure. For Lie groups the existence of such measures is especially simple: Let G be 
a Lie group. Then there is a volume form Ub5, unique up to nonzero multiplicative 
constants, that is left-invariant. If G is compact, Ub5 is right invariant as well. 



4.6 Actions of Lie Groups on Smooth Manifolds 

Let M be a smooth manifold. An action of a Lie group G (with the unit element 
e) on M is a smooth map (j) : G x M ^ M, such that for aW x £ M and g,h £ G, 
(i) (p{e,x) = X and (ii) (p {g, (j){h, x)) = 4){gh,x). In other words, letting (pg : x £ 
M (j)g{x) = (j){g,x) £ M, we have (i') (p^, = idM and (ii') (pg° (ph — 4>gh- 'Pg is a 
diffeomorphism, since (</'g)~^ = We say that the map g £ G ^ cpg £ Dif f{M) 

is a homomorphism of G into the group of diffeomorphisms of Af . In case that M 
is a vector space and each (pg is a linear operator, the function of G on M is called 
a representation of G on M (see, e.g. [H EJ [51 [9] ) . 

An action of G on M is said to be transitive group action, if for every x,y £ M, 
there is g £ G such that (p{g, x) = y; effective group action, if (pg = idM implies g = e, 
that IS g 1-^ (pg is 1-1; and free group action, if for each x £ M, g i— )• (pg{x) is 1-1. 

For example. 
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1. G = M acts on M = M by translations; explicitly, 

(p-.GxM^M, (t){s,x) = X + s. 

Then for x G M, Ox = Hence M/G is a single point, and the action is 
transitive and free. 

2. A complete flow (j)^ of a vector-field X on M gives an action of M on M, namely 

G R X M (j)^{x) G M. 

3. Left translation Lg : G ^ G defines an effective action of G on itself. It is also 
transitive. 

4. The coadjoint action of G on g* is given by 

Ad* : {g,a) e G X Q* ^ A(i*_i(a) = {Te{Rg-i o Lg))* a e g*. 

Let (j) be an action of G on M. For x £ M the orbit of x is defined by 

O. = {(l)g{x)\g eG}cM 
and the isotropy group of at x is given by 

Gx = {g(^ G\(l){g,x) =x} CG. 

An action </> of G on a manifold M defines an equivalence relation on M by the 
relation belonging to the same orbit; explicitly, for x,y G M, we write x ~ y if there 
exists a g £ G such that (p{g, x) = y, that is, ii y G Ox- The set of all orbits M/G is 
called the group orbit space (see, e.g. [H El El [9] ) . 

For example, let M = M^\{0}, G = SO{2), the group of rotations in plane, and 
the action of G on M given by 



(x, y) I I — > (x cos — y sin 9, x sin 6 + y cos ( 



cos t> — sm ti 
sin 9 cos 9 

The action is always free and effective, and the orbits are concentric circles, thus the 
orbit space is M/G ~ R*_^_. 

A crucial concept in mechanics is the infinitesimal description of an action. Let 
(/> : G X M — 7- M be an action of a Lie group G on a smooth manifold M. For each 

C eg, 

(j)^:Rx M ^ M, (j)^{t, x) = (/) (exp(i^), x) 
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is an R — action on M. Therefore, ^exp(t^) : M ^ M is a flow on M; the correspond- 
ing vector-field on M, given by 



^exp(to(2^) 

t=0 



is called the infinitesimal generator of the action, corresponding to in g. 
The tangent space at x to an orbit Ox is given by 

TA = {^Mi^M e 0}- 

Let : G X M — ^ M be a smooth G action. For all g e G, all ^,7] E Q and 

all a, /? G M, we have: 

i^dg^,) ^ = </'*-iCm> [Cm> Vm] = - [C, v]m^ and (a^ + /3rj)M = a^M + I^Vm- 

Let M be a smooth manifold, G a Lie group and (j) : G x M ^ M a G— action 

on M. We say that a smooth map f : M ^ M is with respect to this action if for 

all g eG, 

f °(f>g = (t>g° f- 

Let / : M — )■ M be an equivariant smooth map. Then for any ^ G g we have 

° Cm = Cm o /■ 



4.7 Basic Tables of Lie Groups and Their Lie Algebras 

One classifies Lie groups regarding their algebraic properties (simple, semisimple, 
solvable, nilpotent, Abelian), their connectedness (connected or simply connected) 
and their compactness (see Tables A.1-A.3). This is the content of the Hilbert 5th 
problem. 
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Some real Lie groups and their Lie algebras: 



Lie 

group 


Description 


Remarks 


Lie 
algb. 


Description 


dim 

/-R 


R" 


Euclidean space 
with addition 


Abelian, simply 
connected, not 
compact 


K" 


the Lie bracket 

is zero 


n 


R"" 


nonzero real 
numbers with 

multiplication 


Abelian, not 
connected, not 

compact 


R 


the Lie bracket 
is zero 


1 


E>" 


positive real 
numbers with 
multiplication 


Abelian, simply 
connected, not 
compact 


R 


the Lie bracket 
is zero 


1 


S' = 
R/Z 


complex num- 
bers of absolute 
value 1, with 
multiplication 


Abelian, con- 
nected, not sim- 
ply connected. 
comptUt 


R 


the Lie bracket 
is zero 


1 




noil zero 
quaternions 
with multiplicBr 
tion 


simply con- 
nected, not 
compact 


H 


quaternions, 
with Lie bracket 
the commutator 


4 




quaternions of 
absolute value 
1, with mul- 
tiplication; a 
3— sphere 


simply con- 
nected, com- 
pact, simple 
and semi- 
simple, isomor- 
phic to SU{2), 
SO{3) and to 
Spin{3) 


R^ 


real 3— vectors, 
with Lie bracket 
the cross prod- 
uct; isomorphic 
to 5u(2) and to 
so(3) 


3 


GL{n, R) 


general linear 

group: invert- 
ible n—hy-n 
real matrices 


not connected, 
not compact 


M(n,R; 


n— by-n matri- 
ces, with Lie 
bracket the 
commutator 




GL+(n,]R; 


n—hy-n real 
matrices with 

positive deter- 
minant 


simply con- 
nected, not 
compact 


M(n, E; 


n— by-n matri- 
ces, with Lie 
bracket the 
commutator 


a' 
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Classical real Lie groups and their Lie algebras: 



Lie 

group 


Description 


Remarks 


Lie 
algb. 


Description 


dim 

/-R 


SL{n,R) 


special linear 
group: real 
matrices with 
determinant 1 


simply con- 
nected, not 
compact if 
n> 1 


Sl(n, K) 


square matrices 
with trace 0, 
with Lie bracket 
the commutator 


•1 

n — 
1 


0{n,R) 


orthogonal 
group: real 
orthogonal 
matrices 


not connected, 
compact 


so(n,K; 


skew- 
symmetric 
square real 
matrices, with 
Lie bracket 
the commuta- 
tor; so(3,K) is 
isomorphic to 
5u(2) and to 
with the cross 
product 


n{n- 
l)/2 


SO{n,R) 


special orthogo- 
nal group: real 
orthogonal ma- 
trices with de- 
terminant 1 


connected, com- 
pact, for n > 2: 
not simply con- 
nected, for n — 
3 and n > 
5: simple and 
semisimple 


so(n,E; 


skcw- 
symmetric 
square real 
matrices, with 
Lie bracket the 
commutator 


n{n- 
l)/2 


Spin{n) 


spinor group 


simply con- 
nected, com- 
pact, for n — 3 
and n > 5: 
simple and 
semisimple 


so(n,R; 


skew- 
symmetric 
square real 
matrices, with 
Lie bracket the 
commutator 


n{n- 
l)/2 


Uin) 


unitary group: 
complex unitary 
n—hy-n matri- 
ces 


isomorphic to 
for n = 1, 
not simply 
connected, 
compact 


u(n) 


square com- 
plex matrices 
A satisfying 
A = -A*, with 
Lie bracket the 
commutator 




SU{n) 


special unitary 
group: complex 
unitary n— by-n 
matrices with 
determinant 1 


simply con- 
nected, com- 
pact, for n > 2: 
simple and 
semisimple 


su{n) 


square complex 
matrices A with 

trace satisfy- 
ing A = -A*, 
with Lie bracket 
the commutator 


n'- 
1 
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Basic complex Lie groups and their Lie algebras 




Lie 
group 


Description 


Remarks 


Lie 
algb. 


Description 


dim 

/C 


C" 


group operation 
is addition 


Abelian, simply 
connected, not 
compact 


C" 


the Lie bracket 
is zero 


n 




nonzero com- 
plex numbers 
witli multiplica- 
tion 


Abelian, not 
simply con- 
nected, not 
compact 


C 


the Lie bracket 
is zero 


1 


GL(n, C) 


general linear 
group: invert- 
ible n—hy-n 
complex matri- 
ces 


simply con- 
nected, not 
compact, for 
n — 1: iso- 
morphic to 


M(n,C 


1 n— by-n matri- 
ces, with Lie 
bracket the 
commutator 


n' 


SL{n,C) 


special linear 
group: complex 
matrices with 
determinant 1 


simple, 
semisimple, 
simply con- 
nected, for 
n > 2: not 
compact 


s[(n,C) 


square matrices 
with trace 0, 
with Lie bracket 
the commutator 


n'- 
1 


0{n,C) 


orthogonal 
group: com- 
plex orthogonal 
matrices 


not connected, 
for n > 2: not 
compact 


so(n, 


skew- 
symmetric 
square complex 
matrices, with 
Lie bracket the 
commutator 


n{n 
l)/2 


SO{n,C) 


special orthogo- 
nal group: com- 
plex orthogonal 
matrices with 
determinant 1 


for n > 2: 
not compact, 
not simply 
connected, for 
n = 3 and 
n > 5: simple 
and semisimple 


so(n, 


skew- 
symmetric 
square complex 
matrices, with 
Lie bracket the 
commutator 


n{n 
l)/2 



4.8 Representations of Lie groups 

The idea of a representation of a Lie group plays an important role in the study of 
continuous symmetry (see, e.g., [22]). A great deal is known about such representa- 
tions, a basic tool in their study being the use of the corresponding 'infinitesimal' 
representations of Lie algebras. 

Formally, a representation of a Lie group G on a vector space V (over a field K) 

^'^The dimensions given are dimensions over C. Note that every complex Lie group/algebra can 
also be viewed as a real Lie group/algebra of twice the dimension. 
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is a group homomorphism G — t- AutiV) from G to the automorphism group of V . 
If a basis for the vector space V is chosen, the representation can be expressed as a 
homomorphism into GL{n, K). This is known as a matrix representation. 

On the Lie algebra level, there is a corresponding linear map from the Lie algebra 
of G to Endiy) preserving the Lie bracket [•,•]. 

If the homomorphism is in fact an monomorphism, the representation is said to 
be faithful. 

A unitary representation is defined in the same way, except that G maps to 
unitary matrices; the Lie algebra will then map to skew-Hermitian matrices. 

Now, if G is a semisimple group, its finite-dimensional representations can be 
decomposed as direct sums of irreducible representations. The irreducibles are in- 
dexed by highest weight; the allowable (dominant) highest weights satisfy a suitable 
positivity condition. In particular, there exists a set of fundamental weights, indexed 
by the vertices of the Dynkin diagram of G (see below) , such that dominant weights 
are simply non-negative integer linear combinations of the fundamental weights. 

If G is a commutative compact Lie group, then its irreducible representations 
are simply the continuous characters of G. A quotient representation is a quotient 
module of the group ring. 

4.9 Root Systems and Dynkin Diagrams 

A root system is a special configuration in Euclidean space that has turned out to 
be fundamental in Lie theory as well as in its applications. Also, the classification 
scheme for root systems, by Dynkin diagrams, occurs in parts of mathematics with 
no overt connection to Lie groups (such as singularity theory, see e.g., [22]). 

4.9.1 Definitions 

Formally, a root system is a finite set ^> of non-zero vectors (roots) spanning a 
finite-dimensional Euclidean space V and satisfying the following properties: 

1. The only scalar multiples of a root a in V which belong to $ are a itself and 
-a. 

2. For every root a m V, the set $ is symmetric under refiection through the 
hyperplane of vectors perpendicular to a. 

3. If a and /? are vectors in <I>, the projection of 2/3 onto the line through a is an 
integer multiple of a. 

The rank of a root system $ is the dimension of V. Two root systems may 
be combined by regarding the Euclidean spaces they span as mutually orthogonal 
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subspaces of a common Euclidean space. A root system which does not arise from 
such a combination, such as the systems A2, B2, and G2 in Figure O is said to be 
irreducible. 

Two irreducible root systems {Vi, ^i) and (V2, ^2) are considered to be the same 
if there is an invertible linear transformation Vi ^ V2 which preserves distance up 
to a scale factor and which sends to ^2- 

The group of isometrics of V generated by reflections through hyperplanes as- 
sociated to the roots of $ is called the Weyl group of $ as it acts faithfully on the 
finite set the Weyl group is always finite. 

4.9.2 Classification 

It is not too difficult to classify the root systems of rank 2 (see Figure [3]) . 



AikAi Aj Bs Gs 




Figure 3: Classification of root systems of rank 2. 

Whenever <I> is a root system in V and is a subspace of V spanned by ^' = 
$ n W, then ^ is a root system in W. Thus, our exhaustive list of root systems 
of rank 2 shows the geometric possibilities for any two roots in a root system. In 
particular, two such roots meet at an angle of 0, 30, 45, 60, 90, 120, 135, 150, or 180 
degrees. 

In general, irreducible root systems are specified by a family (indicated by a 
letter A to G) and the rank (indicated by a subscript n). There are four infinite 
families: 

• An{n > 1), which corresponds to the special unitary group, SU{n + 1); 

• Bn (n > 2), which corresponds to the special orthogonal group, S0{2n + 1); 

• C„ (n > 3), which corresponds to the symplectic group, Sp{2n); 

• Dn (n > 4), which corresponds to the special orthogonal group, SO{2n), 

as well as five exceptional cases: Eq,Ey,Es,F4,G2- 
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4.9.3 Dynkin Diagrams 



A Dynkin diagram is a graph with a few different kinds of possible edges (see Figure 
S]). The connected components of the graph correspond to the irreducible subalge- 
bras of 0. So a simple Lie algebra's Dynkin diagram has only one component. The 
rules are restrictive. In fact, there are only certain possibilities for each component, 
corresponding to the classification of semi-simple Lie algebras (see, e.g., [25]). 




Figure 4: The problem of classifying irreducible root systems reduces to the problem 
of classifying connected Dynkin diagrams. 

The roots of a complex Lie algebra form a lattice of rank k in a Cartan subalgebra 
f) C 0, where k is the Lie algebra rank of g. Hence, the root lattice can be considered 
a lattice in M'^. A vertex, or node, in the Dynkin diagram is drawn for each Lie 
algebra simple root, which corresponds to a generator of the root lattice. Between 
two nodes a and /3, an edge is drawn if the simple roots are not perpendicular. One 
line is drawn if the angle between them is 2tt/3, two lines if the angle is 37r/4, and 
three lines are drawn if the angle is 57r/6. There are no other possible angles between 
Lie algebra simple roots. Alternatively, the number of lines between the simple 
roots a and /? is given by 



N = AafsA^a = ' J.. = 4 cos 



-,2 , 



where A^j^ = is an entry in the Cartan matrix {A^is) (for details on Cartan 

matrix see, e.g., |22]). In a Dynkin diagram, an arrow is drawn from the longer root 
to the shorter root (when the angle is 37r/4 or Svr/G). 

Here are some properties of admissible Dynkin diagrams: 

1. A diagram obtained by removing a node from an admissible diagram is admis- 
sible. 

2. An admissible diagram has no loops. 

3. No node has more than three lines attached to it. 
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4. A sequence of nodes with only two single lines can be collapsed to give an 
admissible diagram. 

5. The only connected diagram with a triple line has two nodes. 

A Coxeter-Dynkin diagram, also called a Coxeter graph, is the same as a Dynkin 
diagram, but without the arrows. The Coxeter diagram is sufficient to characterize 
the algebra, as can be seen by enumerating connected diagrams. 

The simplest way to recover a simple Lie algebra from its Dynkin diagram is to 
first reconstruct its Cartan matrix (Aij). The ith node and jth node arc connected 
by AijAji lines. Since A^j = iff Aji = 0, and otherwise A^j G {—3, —2, —1}, it is 
easy to find Aij and Aji, up to order, from their product. The arrow in the diagram 
indicates which is larger. For example, if node 1 and node 2 have two lines between 
them, from node 1 to node 2, then A12 = — 1 and A21 = —2. 

However, it is worth pointing out that each simple Lie algebra can be constructed 
concretely. For instance, the infinite families An, Bn, Cn, and D„ correspond to 
the special linear Lie algebra gl(n + 1, C), the odd orthogonal Lie algebra so(2n + 
1,C), the symplectic Lie algebra sp(2n, C), and the even orthogonal Lie algebra 
so(2n, C). The other simple Lie algebras are called exceptional Lie algebras, and 
have constructions related to the octonions. 

To prove this classification Theorem, one uses the angles between pairs of roots to 
encode the root system in a much simpler combinatorial object, the Dynkin diagram. 
The Dynkin diagrams can then be classified according to the scheme given above. 

To every root system is associated a corresponding Dynkin diagram. Otherwise, 
the Dynkin diagram can be extracted from the root system by choosing a base, that 
is a subset A of $ which is a basis of V with the special property that every vector 
in $ when written in the basis A has either all coefficients > or else all < 0. 

The vertices of the Dynkin diagram correspond to vectors in A. An edge is 
drawn between each non-orthogonal pair of vectors; it is a double edge if they make 
an angle of 135 degrees, and a triple edge if they make an angle of 150 degrees. In 
addition, double and triple edges are marked with an angle sign pointing toward the 
shorter vector. 

Although a given root system has more than one base, the Weyl group acts 
transitively on the set of bases. Therefore, the root system determines the Dynkin 
diagram. Given two root systems with the same Dynkin diagram, we can match up 
roots, starting with the roots in the base, and show that the systems are in fact the 
same. 

Thus the problem of classifying root systems reduces to the problem of classifying 
possible Dynkin diagrams, and the problem of classifying irreducible root systems 
reduces to the problem of classifying connected Dynkin diagrams. Dynkin diagrams 
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encode the inner product on E in terms of the basis A, and the condition that this 
inner product must be positive definite turns out to be all that is needed to get the 
desired classification (see Figured]). 

In detail, the individual root systems can be realized case~by-case, as in the 
following paragraphs: 

An. Let V be the subspace of R"+^ for which the coordinates sum to 0, and 
let ^ be the set of vectors in V of length \/2 and with integer coordinates in M""*"^. 
Such a vector must have all but two coordinates equal to 0, one coordinate equal to 
1, and one equal to -1, so there are + n roots in all. 

B„. Let V = M", and let $ consist of all integer vectors in V of length 1 or \/2. 
The total number of roots is 2n^. 

C^: Let V = M", and let $ consist of all integer vectors in V of together 
with all vectors of the form 2 A, where A is an integer vector of length 1. The total 
number of roots is 2n^. The total number of roots is 2n^. 

Dn. Let V = M", and let $ consist of all integer vectors in V of length \/2. The 
total number of roots is 2n(n — 1). 

E„. For Vg, let 1/ = M^, and let Eg denote the set of vectors a of length \/2 
such that the coordinates of 2a are all integers and are either all even or all odd. 
Then £'7 can be constructed as the intersection of Eg with the hyperplane of vectors 
perpendicular to a fixed root a in Eg, and Eq can be constructed as the intersection 
of Eg with two such hyperplanes corresponding to roots a and (3 which are neither 
orthogonal to one another nor scalar multiples of one another. The root systems 
Eq, Ef, and £^8 have 72, 126, and 240 roots respectively. 

F4. For F4, let V = M.'^, and let $ denote the set of vectors a of length 1 or \/2 
such that the coordinates of 2a are all integers and are either all even or all odd. 
There are 48 roots in this system. 

G2. There are 12 roots in G2, which form the vertices of a hexagram. 

4.9.4 Irreducible Root Systems 

Irreducible root systems classify a number of related objects in Lie theory, notably: 

1. Simple complex Lie algebras; 

2. Simple complex Lie groups; 

3. Simply connected complex Lie groups which are simple modulo centers; and 

4. Simple compact Lie groups. 

In each case, the roots are non-zero weights of the adjoint representation. 

A root system can also be said to describe a plant's root and associated systems. 
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4.10 Simple and Semisimple Lie Groups and Algebras 

A simple Lie group is a Lie group which is also a simple group. These groups, and 
groups closely related to them, include many of the so-called classical groups of 
geometry, which lie behind projective geometry and other geometries derived from 
it by the Erlangen programme of Felix Klein. They also include some exceptional 
groups, that were first discovered by those pursuing the classification of simple Lie 
groups. The exceptional groups account for many special examples and configu- 
rations in other branches of mathematics. In particular the classification of finite 
simple groups depended on a thorough prior knowledge of the 'exceptional' possi- 
bilities. 

The complete listing of the simple Lie groups is the basis for the theory of the 

semisimple Lie groups and reductive groups, and their representation theory. This 
has turned out not only to be a major extension of the theory of compact Lie groups 
(and their representation theory), but to be of basic significance in mathematical 
physics. 

Such groups arc classified using the prior classification of the complex simple 
Lie algebras. It has been shown that a simple Lie group has a simple Lie algebra 
that will occur on the list given there, once it is complexified (that is, made into a 
complex vector space rather than a real one) . This reduces the classification to two 
further matters. 

The groups SO{p, q, M) and SO{p + M), for example, give rise to different real 
Lie algebras, but having the same Dynkin diagram. In general there may be different 
real forms of the same complex Lie algebra. 

Secondly, the Lie algebra only determines uniquely the simply connected (uni- 
versal) cover G* of the component containing the identity of a Lie group G. It may 
well happen that G* is not actually a simple group, for example having a non-trivial 
center. We have therefore to worry about the global topology, by computing the 
fundamental group of G (an Abelian group: a Lie group is an iJ— space). This was 
done by Elic Cart an. 

For an example, take the special orthogonal groups in even dimension. With 
— / a scalar matrix in the center, these are not actually simple groups; and having 
a two-fold spin cover, they aren't simply-connected either. They lie 'between' G* 
and G, in the notation above. 

Recall that a semisimple module is a module in which each submodule is a direct 
summand. In particular, a semisimple representation is completely reducible, i.e., 
is a direct sum of irreducible representations (under a descending chain condition). 
Similarly, one speaks of an Abelian category as being semisimple when every object 
has the corresponding property. Also, a semisimple ring is one that is semisimple as 
a module over itself. 
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A semisimple matrix is diagonalizable over any algebraically closed field contain- 
ing its entries. In practice this means that it has a diagonal matrix as its Jordan 
normal form. 

A Lie algebra q is called semisimple when it is a direct sum of simple Lie algebras, 
i.e., non-trivial Lie algebras £ whose only ideals are {0} and £ itself. An equivalent 
condition is that the Killing form 

B{X, Y) = Ti{Ad{X) Ad{Y)) 

is non-degenerate [27]. The following properties can be proved equivalent for a 
finite-dimensional algebra £ over a field of characteristic 0: 

1. £ is semisimple. 

2. £ has no nonzero Abelian ideal. 

3. £ has zero radical (the radical is the biggest solvable ideal). 

4. Every representation of £ is fully reducible, i.e., is a sum of irreducible repre- 
sentations. 

5. £ is a (finite) direct product of simple Lie algebras (a Lie algebra is called 
simple if it is not Abelian and has no nonzero ideal ) . 

A connected Lie group is called semisimple when its Lie algebra is semisimple; 
and the same holds for algebraic groups. Every finite dimensional representation 
of a semisimple Lie algebra, Lie group, or algebraic group in characteristic is 
semisimple, i.e., completely reducible, but the converse is not true. Moreover, in 
characteristic p > 0, semisimple Lie groups and Lie algebras have finite dimensional 
representations which are not semisimple. An element of a semisimple Lie group or 
Lie algebra is itself semisimple if its image in every finite-dimensional representation 
is semisimple in the sense of matrices. 

Every semisimple Lie algebra g can be classified by its Dynkin diagram |22j . 

5 Some Classical Examples of Lie Groups 
5.1 Galilei Group 

The Galilei group is the group of transformations in space and time that connect 
those Cartesian systems that are termed 'inertial frames' in Newtonian mechanics. 
The most general relationship between two such frames is the following. The origin 
of the time scale in the inertial frame S' may be shifted compared with that in S; 
the orientation of the Cartesian axes in S' may be different from that in S; the 
origin O of the Cartesian frame in S' may be moving relative to the origin O in 5 
at a uniform velocity. The transition from S to S' involves ten parameters; thus the 
Galilei group is a ten parameter group. The basic assumption inherent in Galilei- 
Newtonian relativity is that there is an absolute time scale, so that the only way 
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in which the time variables used by two different 'inertial observers' could possibly 
differ is that the zero of time for one of them may be shifted relative to the zero of 
time for the other (see, e.g. [HEKH]). 

Galilei space-time structure involves the following three elements: 

1. World, as a 4D affine space A^. The points of are called world points or 
events. The parallel transitions of the world A'^ form a linear (i.e., Euclidean) 
space M^. 

2. Time, as a linear map t : — t- M of the linear space of the world parallel 
transitions onto the real 'time axes'. Time interval from the event a G to 
b G is called the number t(b — a); if t{b — a) = then the events a and b are 
called synchronous. The set of all mutually synchronous events consists a 3D 
affine space A\ being a subspace of the world A"^. The kernel of the mapping 
t consists of the parallel transitions of A"^ translating arbitrary (and every) 
event to the synchronous one; it is a linear 3D subspace of the space M^. 

3. Distance (metric) between the synchronous events, 

p{a,b) =\\ a — b \\, for all a,bGA^, 

given by the scalar product in M'^. The distance transforms arbitrary space of 
synchronous events into the well known 3D Euclidean space E^. 

The space A'^, with the Galilei space-time structure on it, is called Galilei space. 
Galilei group is the group of all possible transformations of the Galilei space, pre- 
serving its structure. The elements of the Galilei group are called Galilei transfor- 
mations. Therefore, Galilei transformations are affine transformations of the world 
A^ preserving the time intervals and distances between the synchronous events. 

The direct product M x R^, of the time axes with the 3D linear space R3 with 
a fixed Euclidean structure, has a natural Galilei structure. It is called Galilei 
coordinate system. 

5.2 General Linear Group 

The group of linear isomorphisms of to M" is a Lie group of dimension , called 
the general linear group and denoted G/(n,R). It is a smooth manifold, since it is 
a subset of the vector space L(M",M") of all hnear maps of R" to M", as G/(n,M) 
is the inverse image of R\{0} under the continuous map A i— )• det A of L(M", M") to 
M. The group operation is composition (see, e.g. [U El HI E] ) ■ 

(A, B) £ Gl{n, M)xGl{n,M.)^ AoBe Gl{n, M) 
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and the inverse map is 

A G Gi(n,R) ^ A^^ e Gl{n,R). 

If we choose a basis in R", we can represent each clement A G Gl{n,M) by an 
invertible (n x n) — matrix. The group operation is then matrix multiplication and 
the inversion is matrix inversion. The identity is the identity matrix The group 
operations are smooth since the formulas for the product and inverse of matrices are 
smooth in the matrix components. 

The Lie algebra of G/(n,M) is 0[(n), the vector space L(R"',M") of all linear 
transformations of M", with the commutator bracket 

[A, B]=AB- BA. 

For every ^ G -L(M", M'*), 

°° f 

7^ : i G M ^7^(i) = J] e Gl{n, R) 

is a one-parameter subgroup of GZ(n,R), because 

7^(0) = / and ^,,{t) = Y,j—^,A^ = lA{t)A 

Hence 7^ is an integral curve of the left -invariant vector-field Xa- Therefore, the 
exponential map is given by 



exp : A G L(M",M") ^ exp(^) = 




G G/(n,R). 



1=0 



For each A G GZ(n,R) the corresponding adjoint map 

Ad A ■■ -L(R",R") ^ i:(R",M") 

is given by 

AdAB = A-B-A-'^. 

5.3 Rotational Lie Groups in Human/Humanoid Biomechanics 

Local kinematics at each rotational robot or (synovial) human joint, is defined as 

a group action of an nD constrained rotational Lie group SO{n) on the Euclidean 
space R". In particular, there is an action of SO {2) — group in uniaxial human joints 
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(cylindrical, or hinge joints, like knee and elbow) and an action of 5*0 (3) — group 
in three-axial human joints (spherical, or ball-and-socket joints, like hip, shoulder, 
neck, wrist and ankle). In both cases, SO{n) acts, with its operators of rotation, 
on the vector x = {x'^}, {i = 1, 2, 3) of external, Cartesian coordinates of the parent 
body-segment, depending, at the same time, on the vector q = {g*}, (s = 1, • • • , n) 
on n group-parameters, i.e., joint angles (see O [U [51 [9]). 
Each joint rotation R € SO{n) defines a map 

R-.x^"^ x^, i?(x^, q') = Rgsxf", 

where Rqs £ SO{n) are joint group operators. The vector v = {vg}, (s = 1, • • • ,n) 
of n infinitesimal generators of these rotations, i.e., joint angular velocities, given by 



dR{xf,q'' 



dq^ 



d 



constitute an nD Lie algebra so(n) corresponding to the joint rotation group SO{n). 
Conversely, each joint group operator Rqs , representing a one-parameter subgroup of 
SO{n), is defined as the exponential map of the corresponding joint group generator 

Vs 

Rqs = expiq^'vs) 

This exponential map represents a solution of the joint operator differential equation 
in the joint group-parameter space {q^} 

dRqB 

— = VsRqB ■ 



5.3.1 Uniaxial Group of Joint Rotations 

The uniaxial joint rotation in a single Cartesian plane around a perpendicular axis, 
e.g., xy— plane about the z axis, by an internal joint angle 9, leads to the following 
transformation of the joint coordinates: 

x' = X cos 9 — y sin 9, y' = x sin 9 + y cos 9. 

In this way, the joint S'0(2)— group, given by 

acts in a canonical way on the Euclidean plane by 

Qnfo\ — ( ( ^^^^ —sm9 \ f x \ \ / xcos^ — ysin^ \ 

Y Y ^ cos 9 J \ y J J \ X sin 6 y cos J 
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Its associated Lie algebra 50 (2) is given by 

since the curve 7^ G SO {2} given by 

1 



t G 



COS to — sin t9 
sin to cos 



G 50(2), 



passes through the identity I2 = 



dt 



1 

leit) 



t=o 



and then 

-e 
e 



so that I2 is a basis of so(2), since dim (S'0(2)) = 1. 

The exponential map exp : so (2) SO (2) is given by 



exp 



-9 

e 



= 7e(l) 



cos to —sin to 
sin cos 



The infinitesimal generator of the action of 50(2) on R^, i.e., joint angular 
velocity v, is given by 

d d 
ox oy 



since 



^^R2 {x, y) = ^ 



t=0 



exp(ti;) (x, y) = 



t=o 



cos — sin tti 
sin tv cos ii; 



The momentum map J : T* 



associated to the lifted action of SO (2) on 



by 



J{x,y,Pi,P2) = xpy-ypx, since 
J {x, y,Px,Py) (0 = {Pxdx + pydy){vji^2) = -vp^y + -vpyX. 

The Lie group SO{2) acts on the symplectic manifold (R^, co = dpxAdx+dpyAdx) 



, {x,y,px,Py) 



cos ^ — sin ( 

sin cos 

{x cos ^ — y sin ^, x sin ^ + y cos 0, Px cos — py sin ^, p^; sin +py cos 0) , 
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5.3.2 Three Axial Group of Joint Rotations 



The three-axial 50(3)— group of human-hke joint rotations depends on three pa- 
rameters, Euler joint angles = ip, 9), defining the rotations about the Cartesian 
coordinate triedar (x, y, z) placed at the joint pivot point. Each of the Euler angles 
are defined in the constrained range (— vr, vr), so the joint group space is a constrained 
sphere of radius vr (see El El E] ) . 

Let G = 50(3) = {A e 7W3x3(M) : = /3,det(A) = 1} be the group 

of rotations in W^. It is a Lie group and dim(G) = 3. Let us isolate its one- 
parameter joint subgroups, i.e., consider the three operators of the finite joint rota- 



tions Rip, 



Re G SO{3), given by 



R,, 



1 

cos if — sin If 
sin If cos If 



R,i 



cos ^ sin ip 

1 
— sin ip cos V' 



Rf) 



COS ( 

sin ( 




- sin 6* 
COS 61 
1 



corresponding respectively to rotations about x— axis by an angle ip, about y— axis 
by an angle ip, and about z— axis by an angle 6. 

The total three-axial joint rotation A is defined as the product of above one- 
parameter rotations R^,, Rjj, Rg, i.e., A = Ru, ■ R^ ■ Re is equals 



A 



cos Ip COS ip — cos 6 sin sin ip cos tp cos ip -\- cos 6 cos p sin ip sin 9 sin ip 
— sin ip cos ip — cos 9 sin p sin ip — sin ip sin pi + cos 9 cos ip cos %p sin 9 cos ^ip 
sin 9 sin ip — sin 9 cos ip cos 9 



However, the order of these matrix products matters: different order products give 
different results, as the matrix product is noncommutative product. This is the 
reason why Hamilton's quaternion]^ are today commonly used to parameterize the 
S'0(3)— group, especially in the field of 3D computer graphics. 

The one-parameter rotations R^,R^,Rq define curves in 50(3) starting from 
/I 0\ 

-^3 = 1 . Their derivatives in 99 = 0, V' = and 9 = Q belong to the asso- 

V 1 y 

ciated tangent Lie algebra so (3). That is the corresponding infinitesimal generators 
of joint rotations - joint angular velocities t;^, v^,vg G so(3) - are respectively given 



^^Note that this product is noncommutative, so it really depends on the order of multiplications. 

^^Recall that the set of Hamilton's quaternions H represents an extension of the set of complex 
numbers C. We can compute a rotation about the unit vector, u by an angle 9. The quaternion q 
that computes this rotation is 
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by 



V0 





0-1 

1 

0-10 

1 1 




d d 
oz oy 



d d 
= -X— + y- 



1 

-10 



d d 
ox oz 



dy dx 



Moreover, the elements are linearly independent and so 



50 (3) 





a 
—b 7 



a 



|a,6,7 G 



The Lie algebra so(3) is identified with by associating to each v = {v^,v^,vq) G 



the matrix v G so (3) given by v = 
identities: 

1. u X V = [u, v\; and 

2. u-v = -\Tr:{u-v). 

The exponential map exp : so (3) - 



-a 

a 
—b 7 



Then we have the following 



SO (3) is given by Rodrigues relation 

. .. / ll„ll \ 2 



, , ^ sin k; 1 / sm ■ 

expff ) = / H II II V + - — ip 

\\v\\ 2 \ h 



where the norm ||?;|| is given by 



\\v\\ = vVF+(^^^F+PF- 

The the dual, cotangent Lie algebra so (3)*, includes the three joint angular mo- 
menta p<f,Ptij,Pe G so (3)*. derived from the joint velocities v by multiplying them 
with corresponding moments of inertia. 

Note that the parameterization of 50(3)— rotations is the subject of continuous 
research and development in many theoretical and applied fields of mechanics, such 
as rigid body, structural, and multibody dynamics, robotics, spacecraft attitude 
dynamics, navigation, image processing, etc. 
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5.3.3 The Heavy Top 



Consider a rigid body moving with a fixed point but under the influence of gravity. 
This problem still has a configuration space S0{3), but the symmetry group is only 
the circle group S^, consisting of rotations about the direction of gravity. One says 
that gravity has broken the symmetry from SO{3) to S^. This time, eliminating the 
symmetry mysteriously leads one to the larger Euclidean group SE{3) of rigid 
motion of M^. Conversely, we can start with SE{3) as the configuration space for 
the rigid-body and 'reduce out' translations to arrive at SO{3) as the configuration 
space. The equations of motion for a rigid body with a fixed point in a gravitational 
field give an interesting example of a system that is Hamiltonian. The underlying 
Lie algebra consists of the algebra of infinitesimal Euclidean motions in (see 

The basic phase-space we start with is again T*S0{3), parameterized by Euler 
angles and their conjugate momenta. In these variables, the equations are in canon- 
ical Hamiltonian form. However, the presence of gravity breaks the symmetry, and 
the system is no longer S0{3) invariant, so it cannot be written entirely in terms of 
the body angular momentum p. One also needs to keep track of F, the 'direction of 
gravity' as seen from the body. This is defined by F = A~^k, where k points upward 
and A is the element of 5*0 (3) describing the current configuration of the body. The 
equations of motion are 

where is the body's angular velocity vector, Ii,l2,l3 are the body's principal 
moments of inertia, M is the body's mass, g is the acceleration of gravity, x is the 
body fixed unit vector on the line segment connecting the fixed point with the body's 
center of mass, and I is the length of this segment. 

5.4 Euclidean Groups of Rigid Body Motion 

In this subsection we give description of two most important Lie groups in classical 
mechanics in 2D and 3D, 5^(2) and SE{3), respectively (see [11E119]). 



Pi = -^j^P2P3 + Mgl{T^x^ - 

P2 = %^P3Pi + M^^rV - 

P3 = ^^Yl^PiP2 + Mgl{T\^- 
and r = r X 
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5.4.1 Special Euclidean Group SE{2) in the Plane 

The motion in uniaxial human joints is naturally modelled by the special Euclidean 
group in the plane, SE{2). It consists of all transformations of of the form Az+a, 
where z,a eM?, and 



A G 50(2) = < matrices of the form 



cos U — sin ( 
sin 9 cos 9 



In other words, group SE{2) consists of matrices of the form: 
(Rg, a) = ^ ^ , where a G and Rg is the rotation matrix: 

~ shi^ cos^^ ^ ' ^ 3x3 identity matrix. The inverse (i?6i,a)~^ 

is given by 

CD ^-l ( Re a\ ( R-e -R-ea \ 

=V / J =1 I )■ 

The Lie algebra se(2) of SE{2) consists of 3 x 3 block matrices of the form 



Q Q y, where J=( _^ ^ ), {J = J = -J), 



with the usual commutator bracket. If we identify sc(2) with by the isomorphism 

"^■^ ) e5e(2)^(e,t;)eM3, 
then the expression for the Lie algebra bracket becomes 

[{C,Vl,V2), {C,Wl,W2)] = {0,CV2 ~ CW2,^WI - (vl) = {0,^J^W - (j'^v), 

where V = {vi,V2) and w = {wi,W2)- 
The adjoint group action of 

{Re,a)^1' ;) on (C,^)=("f I 

is given by the group conjugation, 

f Re a \ ( -iJ v\ ( R_e -R-ea \ _ ( -^J iJa + Rev 

lo/jloojlo / /~V0 



46 



or, in coordinates, 

Ad^R^^a){i,v) = {i,iJa + Rev). (2) 

In proving ^ we used the identity RqJ = JRg. Identify the dual algebra, se(2)*, 

/ Ej Q \ 

with matrices of the form ( ^ „ | , via the nondegenerate pairing given by the 

via 



a 







trace of the product. Thus, sc(2)* is isomorphic to 

ese(2)* ^ (//,a) e M^ 



f J 



a 







so that in these coordinates, the pairing between se(2)* and se(2) becomes 

that is, the usual dot product in M^. The coadjoint group action is thus given by 

a)~^ ^) ~ if^ ~ ^sot ■ J a + Reo). (3) 

Formula ([3]) shows that the coadjoint orbits are the cylinders r*5^ = {(/x, a)| ||a|| = 
const} if a 7^ together with the points are on the //—axis. The canonical cotangent 
bundle projection vr : T*^^ — )• 5^ is defined as 7r(/i,a) = a. 

5.4.2 Special Euclidean Group 5^(3) in the 3D Space 

The most common group structure in human biodynamics is the special Euclidean 
group in 3D space, SE{3). Briefly, the Euclidean 5£'(3)— group is defined as a 
semidirect (noncommutative) product of 3D rotations and 3D translations, SE(3) := 
S0{3) \> (see [HElin]). Its most important subgroups are the following: 



Subgroup 



Definition 



S0{3), group of rotations 
in 3D (a spherical joint) 



Set of all proper orthogonal 
3 X 3 — rotational matrices 



Set of all 3 X 3 — matrices: 



SE{2), special Euclidean group 
in 2D (all planar motions) 



cos t> 


smU 


rx 


— sin 9 


cos 9 










1 



SO{2), group of rotations in 2D 
subgroup of 5£^(2)-group 
(a revolute joint) 



W, group of translations in 3D 
(all spatial displacements) 



Set of all proper orthogonal 
2 X 2 — rotational matrices 
included in SE{2) — group 



Euclidean 3D vector space 
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Lie Group SE{3) and Its Lie Algebra An element of SE{3) is a pair {A, a) 
where A £ 50(3) and a G M^. The action of SE{3) on is the rotation A followed 
by translation by the vector a and has the expression 



Using this formula, one sees that multiplication and inversion in SE{3) are given by 

{A, a){B, b) = {AB, Ah + a) and {A, a)~^ = {A-^ , -A'^a), 

for A,B £ 50(3) and a, 6 € M?. The identity element is (/, 0). 

The Lie algebra of the Euclidean group SE{3) is se(3) = M'^ x M'^ with the Lie 
bracket 



The Lie algebra of the Euclidean group has a structure that is a special case of 
what is called a semidirect product. Here it is the product of the group of rotations 
with the corresponding group of translations. It turns out that semidirect products 
occur under rather general circumstances when the symmetry in T*G is broken (see 



The dual Lie algebra of the Euclidean group SE{3) is 5e(3)* = x with the 
same Lie bracket ([5]). 

Representation of SE{3) In other words, SE{3) := 50(3) oM'^ is the Lie group 
consisting of isometrics of M^. 

Using homogeneous coordinates, we can represent SE{3) as follows. 



with the action on given by the usual matrix-vector product when we identify 
with the section M'^ x {1} C M^. In particular, given 



{A, a) ■ X = Ax + a. 



[{C, u), {r],v)] = X r],^ X V - rj X u). 



(4) 



laiHii)- 





and (; G M^, we have 



9 ■ q 



Rq+p, 



or as a matrix-vector product 




48 



Lie algebra of SE{3) The Lie algebra of SE{3) is given by 



se(3)= Q ^ j eM^{R):ueso{3),ve 
where the attitude matrix u; : ^ so(3) is given by 





( « 

























The exponential map of SE{?>) The exponential map, exp : se(3) SE{3), is 
given by 

u V \ f expfcj) Av 
exp ' ' — ' 

where 



y V 1 



1 — cos||a;|| llwll — sin IIcjII o 

^ = ^ + — r^^^+ II ii3 ^ ' 



and exp(a;) is given by the Rodriguez' formula, 

sin ||a;|| 1 — cos ||a;|| o 
exp(C(;) = i H t — ^u; H 5 ■ 

ll^ll ||<^|| 

In other words, the special Euclidean group SE{3) := S0{3) \> is the Lie 
group consisting of isometrics of the Euclidean 3D space M'^. An element of SE{3) 
is a pair {A, a) where A G SO{3) and a G M^. The action of SE{3) on is the 
rotation A followed by translation by the vector a and has the expression 

{A, a) ■ X = Ax + a. 

The Lie algebra of the Euclidean group SE{3) is se(3) = x M'^ with the Lie 
bracket 

[i^,y'),{v,v)] = {^xv,^xv-vxu). (5) 
Using homogeneous coordinates, we can represent SE{3) as follows, 

SEi3)= ^ j gGL(4,M) :i?G50(3),pG]B 

with the action on given by the usual matrix-vector product when we identify 
with the section x {1} C M^. In particular, given 
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and g G M^, we have 



or as a matrix-vector product, 



g-q = Rq + p, 



R p \ f q 
1 j I 1 



Rq+p 
1 



The Lie algebra of SE{3), denoted se(3), is given by 



se(3) 



UJ V 





e M4(M) : io G so(3), V G 



where the attitude (or, angular velocity) matrix oj iM.^ ^ so (3) is given by 









Uy 





















The exponential map, exp : se(3) — >■ S'£^(3), is given by 
exp 



where 



A = I + 



Ijj V 





1 — cos IIcjI 



exp(a;) Av 
1 



-UJ + 



\uj — sm a; I 



-UJ 



\\UJ\\ \\L0\\ 

and exp(a;) is given by the Rodriguez' formula, 

sin \\uj\\ 1 — cos 



exp(a;) = I -\ — - — ui + 



-to 



5.5 Basic Mechanical Examples 
5.5.1 S'£;(2) -Hovercraft 

Configuration manifold is {0,x,y) G SE{2), given by matrix 





COS 9 


sin 6 


X 


p = 


— sinO 


cos 9 


y 










1 
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Kinematic equations of motion in Lie algebra st(2): 



P = P 



Q 00 Vx 

— OJ Vy 





= 9,Vx = x,Vy = y). 



Kinetic energy: 



where m, I arc mass and inertia moment of the hovercraft. 
Dynamical equations of motion: 



mvx 
mvy 
Ilo 



mojVy + ui, 
—mcovx + U2, 

TU2, 



where r = —h is the torque applied at distance h from the center-of-mass, while 
ui,U2 are control inputs. 

5.5.2 SO (3) -Satellite 

Configuration manifold is rotation matrix R G SO (3), with associated angular- 
velocity (attitude) matrix u; = (a;i,a;2,'*^3) € 50 (3) « given by 



a; G so(3) 



— ^3 0J2 



ws 

-0O2 Ull 



-UJi 





Kinematic equation of motion in so (3): 

R = Ru, 

Kinetic energy: 

1 rj, 

Ek = 2'*' 

where inertia tensor I is given by diagonal matrix, 

I = diag{Ii, 12,13} ■ 
Dynamical Euler equations of motion: 

Id; = la; X u> + r^u*, 

where x is the cross-product in 3D, Tj are three external torques and = n'(t) are 
control inputs. 
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5.5.3 S'£;(3)-Submarine 

The motion of a rigid body in incompressible, irrotational and inviscid fluid is de- 
fined by the configuration manifold SE{3), given by a pair of rotation matrix and 
translation vector, {R,p) € SE{3), such that angular velocity (attitude) matrix and 
linear velocity vector, (w, v) G 5e(3) « M^. 

Kinematic equations of motion in se(3): 

p = Rv, R = Roj. 

Kinetic energy (symmetrical): 

Ek = ^v^Mv + ^w^Iw, 

where mass and inertia matrices arc diagonal (for a neutrally buoyant ellipsoidal 
body with uniformly distributed mass), 

M = diag{mi,m2,m^}, 
I = diag{hj2,h}- 

Dynamical Kirchhoff equations of motion read: 

Mv = Mv X uj, luJ = la; X a; + Mv x v. 

By including the body-fixed external forces and torques, fi,Ti, with input controls 
= u^t), the dynamical equations become: 

Mv = Mv X oj + fiu\ 
loj = la; X a; + Mv X V + TjU*. 

5.6 Newton— Euler £'(3)— Dynamics 

5.6.1 SO{3) : Euler Equations of Rigid Rotations 

Unforced Euler equations read in vector form 

TT = lo; = TT X a;, with I = diag{Ii, 12,13} 
and in scalar form 

huji = {h - ^3)^2^3 
I2CJ2 = ih - /l)w3Wl. 

/3a;3 = (/i - l2)uiU2 
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Using rotational kinetic-energy Lagrangian 

L{oj) = El"* = h^^luj = ^{hi^l + hi^l + hojl) (*= 'transpose') 

Regarding the angular momentum tt = d^L = lu; = (/iwi, /2CJ2, -^3^3) as a vector, 
we can derive unforced Euler equations: tt = tt x a; as a system of Euler-Lagmnge- 
Kirchhoff equations 

-^d^L = d^L X uj. 
Forced Euler equations read in vector form 

TT + W X TT = T 

and in scalar form 

IiLOi + (/a - /2)w2a;3 = Ti 
/2W2 + {h - h)^2,i^i = T2 
/awa + {I2 - Ii)oJiOJ2 = T3 

5.6.2 SE{3) : Coupled Newton Euler Equations 

Forced coupled Newton-Euler equations read in vector form 

p = Mv = F + p X a;, with M. = diag{mi,m2,ms} 
7r = Ia; = T + 7rxa; + pxv, 1 = diag{Ii, 12,13}, 

with principal inertia moments given in Cartesian coordinates {x, y, z) by density 
p— dependent volume integrals 



h = 



p{z^+y^)dxdydz, h = JJJ p{x^ +y^)dxdydz , I3 = JJJ p{x^ +y^)dxdydz , 



In tensor form, the forced-coupled Newton-Euler equations read 
where the permutation symbol efu is defined as 



+1 if is (1,2, 3), (3, 1,2) or (2,3,1), 

-1 if (i,j, fe) is (3, 2,1), (1,3, 2) or (2,1,3), 
otherwise: i = j 01 j = k or k = i. 



53 



In scalar form these equations read 



P2 = F2 + m^vzoji — miviu!3 
Pg = F^ — m2V20Ji + miviu!2 

TTj = Ti + (m2 - ms)v2Vs + {h - 13)^2^3 
TTa = ^2 + (ma - mi)i;ii;3 + (/a - /i)a;ia;3 
TTg = Ta + (mi — m2)viV2 + {h — I2)ooiuj2 

These coupled rigid-body equations can be derived from the Newton-Euler kinetic 
energy 

Eu = -v*Mv + -aj*Iu; 
"2 2 

or, in tensor form 

E = ^M^ji'v^ + hijcj'cj^. 
Using the Kirchhoff-Lagrangian equations 

^d^Ek = d^EkXu + F 
dt 

-jidwEk = d^Ek X OJ + d^Ek x v + T, 
dt 



or, in tensor form 



-d,.E = 4,id,jE)u'' + Fi, 



d_ 
dt 

j^d^^E = 4,{d,,E)u^'^ + el,{d,^E)v^ + Ti 



we can derive linear and angular momentum covectors 



or, in tensor form 
and in scalar form 



Pi = d^iE, TTi = d^iE, 



P = \Pl,P2,P3] = ['miVi,m2V2,m2V3] 

TT = [7ri,7r2,7ra] = [Iiu;i,/2a;2,/3W3], 



54 



with their respective time derivatives, in vector form 



or, in tensor form 

dt^^ dt ^ ' ^ dt ^ dt 

and in scalar form 



Pi = —Pi = —Oy^E, TTi = — TTi = —o^iE, 



P = \Pl,P2,P3] = ['miVl,m2V2,m3V3] 
TT = [71-1,71-2,71-3] = [/ia;i,/2W2,/3^3]- 

In addition, for the purpose of biomechanical injury prediction/prevention, we 
have hnear and angular jolts, respectively given in vector form by 

F = p — pxcj — pxcj 

T = TT — TT XCJ — TTXcj — pXV — pXV, 

or, in tensor fornj^ 

Fi = Pi-elkPju;'' -el^pjilj'', 
Ti = ^i -ei^j^"" -4k^j^^ -e^PiV^ -ey,v\ 
where the linear and angular jolt covectors are 

F = Fi = M.w = Miji)^ = [Fi,F2,F3], 

T = fi = Iu>=Ii,U^ = [fi,f2,f3], 

where 

are linear and angular jerk vectors. 

In scalar form, the S'i?(3)— jolt expands as 

= Pl - m2U}3V2 + {0J2V3 + ^3^2) - 771,2^2^3, 

= p2 + miD^vi - msojivs - m^vsuji + miViUJs, 
= - miU}2Vi + m20JiV2 - V2ioi - mividj2, 

-- Tfi - (m2 - ma) {V3V2 + V2V3) - {h - h) {'■^?.'^2 + W2W3) , 
= 7i-2 + (mi - ms) {v^vi + viv^) + {h - h) (w3tji + W1W3) , 
= 71-3 - (mi - m2) {v2Vi + ^1^2) - (-^1 - h) (^2(^1 + OJ1CJ2) ■ 



In this paragraph the overdots actuaUy denote the absolute Bianchi (covariant) derivatives, so 
that the jolts retain the proper covector character, which would be lost if ordinary time derivatives 
are used. However, for simplicity, we stick to the same notation. 
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5.7 Symplectic Group in Hamiltonian Mechanics 



Here we give a brief description of symplectic group (see [H El E] ) . 

^ ' , with / the n X n identity matrix. Now, A E 



Let J 



-/ 



is cahed a symplectic matrix if A^J A = J. Let Sp{2n,M.) be the set of 2n x 2n 
symplectic matrices. Taking determinants of the condition A^J A = J gives det A = 
±1, and so A G GL(2n,R). Furthermore, iiA,B€ Sp{2n,R), then {AB)'^J{AB) = 
B'^A^JAB = J. Hence, AB E 5p(2n,R), and \iA^JA = J, then J A = {A^y^J = 
(A-i)^J, so J = (A - 1)^ J^-i, or G Sp{2n,R). Thus, S'p(2n,M) is a group. 
The symplectic Lie group 

Sp{2n,R) = {AG GL(2n,R) : A^JA = j} 
is a noncompact, connected Lie group of dimension 2n^ + n. Its Lie algebra 
sp(2n,M) = {A£ L(M2",M2n^) . ^Tj^ ^ j ^ q| ^ 



called the symplectic Lie algebra, consists of the 2n x 2n matrices A satisfying 
A^JA = 0. 

Consider a particle of mass m moving in a potential V{q), where = {q^, q'^, q"^) G 
M^. Newtonian second law states that the particle moves along a curve q{t) in in 
such a way that mq^ = — gradV^(g*). Introduce the 3D-momentum pi = mq'^, and 
the energy (Hamiltonian) 



1 ^ 



Then 



dH 

dq^ 
dH 
dpi 

and hence Newtonian law 



Now, writing z = {q^,Pi), 
J gra,dH{z) - 



dV 
dq^ 
1 

—Pi 
m 



i=l 



-mq 



-Pi 



and 



F 



mq 



1,2,3), 



is equivalent to Hamiltonian equations 



dH 
dpi' 



Pi 



dH 
dq^ ' 




{q\pi) 
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so the complex Hamiltonian equations read 

z = J gradH{z). 

Now let / : X — > X and write w = f{z). If z{t) satisfies the complex 
Hamiltonian equations then w{t) = f{z{t)) satisfies w = A'^z, where = [dw"^ /dz^] 
is the Jacobian matrix of /. By the chain rule, 

■w = A^JgradH{z) = J A gr ad H{z{w)). 

z w 

Thus, the equations for w{t) have the form of Hamiltonian equations with energy 
K{w) = H{z{w)) iff A^J A = J, that is, iff A is symplectic. A nonlinear transfor- 
mation / is canonical iff its Jacobian matrix is symplectic. Sp{2n, M) is the linear 
invar iance group of classical mechanics. 

6 Medical Applications: Prediction of Injuries 

6.1 General Theory of Musculo— Skeletal Injury Mechanics 

The prediction and prevention of traumatic brain injury, spinal injury and musculo- 
skeletal injury is a very important aspect of preventive medical science. Recently, 
in a series of papers [28l [29l [30], we have proposed a new coupled loading-rate 
hypothesis as a unique cause of all above injuries. This new hypothesis states that 
the unique cause of brain, spinal and musculo-skeletal injuries is a Euclidean Jolt, 
which is an impulsive loading that strikes any part of the human body (head, spine 
or any bone/joint) - in several coupled degrees-of- freedom simultaneously. It never 
goes in a single direction only. Also, it is never a static force. It is always an 
impulsive translational and/or rotational force coupled to some mass eccentricity. 
This is, in a nutshell, our universal Jolt theory of all mechanical injuries. 

To show this, based on the previously defined covariant force law, we have firstly 
formulated the fully coupled Newton-Euler dynamics of: 

1. Brain's micro- motions within the cerebrospinal fluid inside the cranial cavity; 

2. Any local inter-vertebral motions along the spine; and 

3. Any local joint motions in the human musculo-skeletal system. 

Then, from it, we have defined the essential concept of Euclidean Jolt, which 
is the main cause of all mechanical injuries. The Euclidean Jolt has two main 
components: 

1. Sudden motion, caused either by an accidental impact or slightly distorted 
human movement; and 
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Figure 5: Human brain and its S* £'(3)— group of microscopic three-dimensional mo- 
tions within the cerebrospinal fluid inside the cranial cavity. 



2. Unnatural mass distribution of the human body (possibly with some added 
masses), which causes some mass eccentricity from the natural physiological body 
state. 

What does this all mean? Wc will try to explain it in "plain English". As we 
live in a 3D space, one could think that motion of any part of the human body, ei- 
ther caused by an accidental impact or by voluntary human movement, "just obeys 
classical mechanics in 6 degrees-of-freedom: three translations and three rotations" . 
However, these 6 degrees-of-frcedom are not independent motions as it is suggested 
by the standard term "degrees-of-freedom". In reality, these six motions of any 
body in space are coupled. Firstly, three rotations are coupled in the so-called ro- 
tation group (or matrix, or quaternion). Secondly, three translations are coupled 
with the rotation group to give the full Euclidean group of rigid body motions in 
space. A simple way to see this is to observe someone throwing an object in the 
air or hitting a tennis ball: how far and where it will fly depends not only on the 
standard "projectile" mechanics, but also on its local "spin" around all three axes 
simultaneously. Every golf and tennis player knows this simple fact. Once the spin 
is properly defined we have a "fully coupled Newton-Euler dynamics" - to start with. 
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Figure 6: Human body representation in terms of SE(3)/SE(2)-groups of rigid-body 
motion, with the vertebral column represented as a chain of 26 flexibly-coupled 
SE(3)-groups. 

The covariant force law for any biodynamical system (which we introduced earlier 
in our biodynamics books and papers, see our references in the cited papers above) 
goes one step beyond the Newton-Euler dynamics. It states: 

Euclidean Force covector field = 
Body mass distribution x Euclidean Acceleration vector field 

This is a nontrivial biomechanical generalization of the fundamental Newton's 
definition of the force acting on a single particle. Unlike classical engineering me- 
chanics of multi-body systems, this fundamental law of biomechanics proposes that 
forces acting on a multi-body system and causing its motions are fundamentally dif- 
ferent physical quantities from the resulting accelerations. In simple words, forces 
are massive quantities while accelerations are massless quantities. More precisely, 
the acceleration vector field includes all linear and angular accelerations of individual 
body segments. When we couple them all with the total body's mass-distribution 
matrix of all body segments (including all masses and inertia moments), we get the 
force co-vector field, comprising all the forces and torques acting on the individual 
body segments. In this way, we have defined the 6-dimensional Euclidean force for 
an arbitrary biomechanical system. 

Now, for prediction of injuries, we need to take the rate-of-change (or derivative. 
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Figure 7: Schematic latero-frontal view of the left knee joint. Although designed 
to perform mainly flexion/extension (strictly in the sagittal plane) with some re- 
stricted medial/lateral rotation in the semi-flexed position, it is clear that the knee 
joint really has at least six-degrees-of-freedom, including three micro-translations. 
The injury actually occurs when some of these microscopic translations become 
macroscopic, which normally happens only after an external jolt. 



with respect to time) of the Euclidean biomechanical force defined above. In this way, 
we get the Euclidean Jolt, which is the sudden change (in time) of the 6-dimensional 
Euclidean force: 

Euclidean Jolt covector field = 
Body mass distribution x Euclidean Jerk vector field 

And again, it consists of two components: (i) massless linear and angular jerks 
(of all included body segments), and (ii) their mass distribution. For the sake of sim- 
plicity, we can say that the mass distribution matrix includes all involved segmental 
masses and inertia moments, as well as "eccentricities" or "pathological leverages" 
from the normal physiological state. 

Therefore, the unique cause of all brain, spine and musculo-skeletal injuries has 
two components: 

1. Coupled linear and angular jerks; and 

2. Mass distribution with "eccentricities" . 
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In other words, there are no injuries in static conditions without any mass 
eccentricities; all injuries are caused by mutually coupled linear and an- 
gular jerks, which are also coupled with the involved human mass distri- 
bution. 

The Euchdean Jolt causes two forms of discontinuous brain, spine or musculo- 
skeletal injury: 

1. Mild rotational disclinations; and 

2. Severe translational dislocations (or, fractures). 

In the cited papers above, we have developed the soft-body dynamics of biome- 
chanical disclinations and dislocations, caused by the Euclidean Jolt, using the 
Cosserat multipolar viscoelastic continuum model. 

Implications of the new universal theory are various, as follows. 

A. The research in traumatic brain injury (TBI, see Figure [5]) has so far iden- 
tified the rotation of the brain-stem as the main cause of the TBI due to various 
crashes/impacts. The contribution of our universal Jolt theory to the TBI research 
is the following: 

1. Rigorously defined this brain rotation as a mechanical disclination of the 
brain-stem tissue modelled by the Cosserat multipolar soft-body model; 

2. Showing that brain rotation is never uni-axial but always three-axial; 

3. Showing that brain rotation is always coupled with translational disloca- 
tions. This is a straightforward consequence of our universal Jolt theory. 

These apparently 'obvious' facts are actually radically new: we cannot separately 
analyze rapid brain's rotations from translations, because they are in reality always 
coupled. 

One practical application of the brain Jolt theory is in design of helmets. Briefly, 
a 'hard' helmet saves the skull but not the brain; alternatively, a 'soft' helmet pro- 
tects the brain from the collision jolt but does not protect the skull. A good helmet 
is both 'hard' and 'soft'. A proper helmet would need to have both a hard external 
shell (to protect the skull) and a soft internal part (that will dissipate the energy 
from the collision jolt by its own destruction, in the same way as a car saves its 
passengers from the collision jolt by its own destruction). 

Similarly, in designing safer car air-bags, the two critical points will be (i) their 
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placement within the car, and (ii) their "soft-hard characteristics", similar to the 
helmet characteristics described above. 

B. In case of spinal injury (see Figm'e[6]), the contribution of our universal Jolt 
theory is the following: 

1. The spinal injury is always localized at the certain vertebral or inter- vertebral 
point; 

2. In case of severe translational injuries (vertebral fractures or discus herniae) 
they can be identified using X-ray or other medical imaging scans; in case of micro- 
scopic rotational injuries (causing the back-pain syndrome) they cannot be identified 
using current medical imaging scans; 

3. There is no spinal injury without one of the following two causes: 

a. Impulsive rotational -|- translational loading caused by either fast human 
movements or various crashes/impacts; and/or 

b. Static eccentricity from the normal physiological spinal form, caused by 
external loading; 

c. Any spinal injury is caused by a combination of the two points above: 
impulsive rotational -|- translational loading and static eccentricity. 

This is a straightforward consequence of our universal Jolt theory. We cannot 
separately analyze translational and rotational spinal injuries. Also, there are no 
"static injuries" without eccentricity. Indian women have for centuries carried bulky 
loads on their heads without any spinal injuries; they just prevented any load ec- 
centricities and any jerks in their motion. 

The currently used "Principal loading hypothesis" that describes spinal injuries 
in terms of spinal tension, compression, bending, and shear, covers only a small 
subset of all spinal injuries covered by our universal Jolt theory. To prevent spinal 
injuries we need to develop spinal jolt awareness: ability to control all possible im- 
pulsive spinal loadings as well as static eccentricities. 

C. In case of general musculo-skeletal injury (see Figure [7] for the particular case 
of knee injury), the contribution of our universal Jolt theory is the following: 

1. The injury is always localized at the certain joint or bone and caused by an 
impulsive loading, which hits this particular joint/bone in several coupled degrees- 
of-freedom simultaneously; 

2. Injury happens when most of the body mass is hanging on that joint; for 
example, in case of a knee injury, when most of the body mass is on one leg with 
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a semi-flexed knee — and then, caused by some external shock, the knee suddenly 
"jerks" (this can happen in running, skiing, and ball games, as well as various 
crashes/impacts); or, in case of shoulder injury, when most of the body mass is 
hanging on one arm and then it suddenly jerks. 

To prevent these injuries we need to develop musculo-skeletal jolt awareness. For 
example, never overload a flexed knee and avoid any kind of uncontrolled motions 
(like slipping) or collisions with external objects. 

6.2 Analytical Mechanics of Traumatic Brain Injury (TBI) 
6.2.1 The 5^(3)-joh: the cause of TBI 

In this subsection we give a brief on TBI mechanics. For more details and references, 
see [28]. 

In the language of modern dynamics, the microscopic motion of human brain 
within the skull is governed by the Euclidean SE(3)-group of 3D motions (see next 
subsection). Within brain's SE(3)-group we have both SE(3)-kinematics (consisting 
of SE(3)-velocity and its two time derivatives: SE(3)-acceleration and SE(3)-jerk) 
and SE(3)-dynamics (consisting of SE(3)-momentum and its two time derivatives: 
SE(3)-force and SE(3)-jolt), which is brain's kinematics x brain's mass-inertia 
distribution. 

Informally, the external SE(3)-joliP^ is a sharp and sudden change in the SE(3)- 
force acting on brain's mass-inertia distribution (given by brain's mass and inertia 
matrices). That is, a 'delta'-change in a 3D force-vector coupled to a 3D torque- 
vector, striking the head-shell with the brain immersed into the cerebrospinal fluid. 
In other words, the SE(3)-jolt is a sudden, sharp and discontinues shock in all 
6 coupled dimensions of brain's continuous micro-motion within the cerebrospinal 
fluid (Figure [5]), namely within the three Cartesian (x, y, 2:)-translations and the 
three corresponding Euler angles around the Cartesian axes: roll, pitch and yaw. 
If the SE(3)-jolt produces a mild shock to the brain (e.g., strong head shake), it 
causes mild TBI, with temporary disabled associated sensory- motor and/or cognitive 
functions and affecting respiration and movement. If the SE(3)-jolt produces a hard 

^''The mechanical SE(3)-jolt concept is based on the mathematical concept of higher-order tan- 
gency (rigorously defined in terms of jet bundles of the head's configuration manifold), as follows: 
When something hits the human head, or the head hits some external body, we have a collision. 
This is naturally described by the SE(3)"momentum, which is a nonlinear coupling of 3 linear New- 
tonian momenta with 3 angular Eulerian momenta. The tangent to the SE(3)-momentum, defined 
by the (absolute) time derivative, is the SE(3)-force. The second-order tangency is given by the 
SE(3)-jolt, which is the tangent to the SE(3)-force, also defined by the time derivative. 
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shock (hitting the head with external mass), it causes severe TBI, with the total 
loss of gesture, speech and movement. 

The SE(3)-jolt is the absolute time-derivative of the covariant force 1-form (or, 
CO- vector field). The fundamental law of biomechanics is the covariant force law: 

Force co-vector field = Mass distribution x Acceleration vector-field, 

which is formally written (using the Einstein summation convention, with indices 
labelling the three Cartesian translations and the three corresponding Euler angles): 

Ffi = mf,^a'', = 1,...,6) 

where denotes the 6 covariant components of the external "pushing" SE(3)-force 
CO- vector field, rUfn, represents the 6x6 covariant components of brain's inertia- 
metric tensor, while a'^ corresponds to the 6 contravariant components of brain's 
internal SE(3)-acceleration vector-field. 

Now, the covariant (absolute, Bianchi) time-derivative ^(•) of the covariant 
SE(3)-force defines the corresponding external "striking" SE(3)-jolt co-vector 
field: 

£(F^) = m^.j^{an = m,. (a^ + T^^^a'^a>) , (6) 

where %{oy) denotes the 6 contravariant components of brain's internal SE(3)-jerk 
vector-field and overdot (') denotes the time derivative. FJ^^ are the Christoffel's 
symbols of the Levi-Civita connection for the SE(3)-group, which are zero in case 
of pure Cartesian translations and nonzero in case of rotations as well as in the 
full-coupling of translations and rotations. 

In the following, we elaborate on the SE(3)-jolt concept (using vector and tensor 
methods) and its biophysical TBI consequences in the form of brain's dislocations 
and disclinations. 



6.2.2 S'£'(3)— group of brain's micro— motions within the CSF 

The brain and the CSF together exhibit periodic microscopic translational and ro- 
tational motion in a pulsatile fashion to and from the cranial cavity, in the fre- 
quency range of normal heart rate (with associated periodic squeezing of brain's 
ventricles). This micro-motion is mathematically defined by the Euclidean (gauge) 
S'^;(3) -group. 

In other words, the gauge £'(3)— group of Euclidean micro-motions of the brain 
immersed in the cerebrospinal fluid within the cranial cavity, contains matrices of 

the form ^ ^ ) ' ^'^^^^ ^ brain's 3D micro-translation vector and R is brain's 

3D rotation matrix, given by the product R = R,^ ■ ■ Rg of brain's three Eulerian 
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micro-rotations, roll = R^, pitch = R^, yaw = Rg, performed respectively about 
the X— axis by an angle (p, about the y— axis by an angle ■0, and about the z— axis 
by an angle 6, 



1 

cos Lp — sin 
sin ip cos tp 



cos -0 sin -0 

1 
— sin cos '0 



, Re^ 



cos ( 
sint 




-sine 
cos 6* 
1 



Therefore, brain's natural 5* £'(3)— dynamics within the cerebrospinal fluid is 
given by the coupling of Newtonian (translational) and Eulerian (rotational) equa- 
tions of micro-motion. 



6.2.3 Brain's natural S'ii^ (3)— dynamics 

To support our coupled loading-rate hypothesis, we formulate the coupled Newton- 
Euler dynamics of brain's micro-motions within the scull's S'S(3)— group of motions. 
The forced Newton-Euler equations read in vector (boldface) form 



Newton : p = Mv = F -|- p x cj, 

Euler : 7r = Icj = T + 7rxa;-|-pxv, 



(7) 



where x denotes the vector cross product 1^ 

M = Mij = diag{mi,m2,m3} and 1 = lij = diag{Ii,l2,l3}, (i,j = 1,2,3) 

are brain's (diagonal) mass and inertia matrices^ defining brain's mass-inertia 
distribution, with principal inertia moments given in Cartesian coordinates {x, y, z) 
by volume integrals 



h 



p{z^ +y^)dxdydz , I2 



p{x'^+z'^)dxdydz, I3 



p{x^+y'^)dxdydz, 



^^Recall that the cross product u x v of two vectors u and v equals u x v = uvsinOn, where 9 
is the angle between u and v, while n is a unit vector perpendicular to the plane of u and v such 
that u and v form a right-handed system. 

^^In reality, mass and inertia matrices (M, I) are not diagonal but rather full 3x3 positive- 
definite symmetric matrices with coupled mass- and inertia-products. Even more realistic, fully- 
coupled mass-inertial properties of a brain immersed in (incompressible, irrotational and inviscid) 
cerebrospinal fluid are defined by the single non-diagonal 6x6 positive-definite symmetric mass- 
inertia matrix Mse(3), the so-called material metric tensor of the 5'i5(3)— group, which has all 
nonzero mass-inertia coupling products. In other words, the 6x6 matrix A^s_b(3) contains; (i) 
brain's own mass plus the added mass matrix associated with the fluid, (ii) brain's own inertia plus 
the added inertia matrix associated with the potential flow of the fluid, and (iii) all the coupling 
terms between linear and angular momenta. However, for simplicity, in this paper we shall consider 
only the simple case of two separate diagonal 3x3 matrices (M, I). 
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dependent on brain's density p = p{x, y, z), 
V = = [wi, ^2, fs]* and 



U)' = [uJi,UJ2,UJ3f 



(where \ Y denotes the vector transpose) are brain's hnear and angular velocity 
vectora^^l (that is, column vectors), 



and 



T, = [Ti,T2,T3 



are gravitational and other external force and torque co- vectors (that is, row vectors) 
acting on the brain within the scull, 



p = Pi = Mv = [pi,P2,P3] = [miVi,m2V2,m2V2] 

77 = TTi = Iu; = [tTi,-K2, its] = [IiUJi,l2^^2,h^3] 

are brain's linear and angular momentum co-vectors. 

In tensor form, the forced Newton-Euler equations ([7]) read 



and 



Pi 



hju^ =Ti + eli^TTjU + el^pjv 



{i,j,k = 1,2,3) 

k 



where the permutation symbol e^,^ is defined as 



+1 iiii,j,k) is (1,2, 3), (3, 1,2) or (2,3,1), 
-1 if (^,j,A;) is (3, 2,1), (1,3, 2) or (2,1,3), 
otherwise: i = j or j = k or k = i. 



In scalar form, the forced Newton-Euler equations ([7D expand as 

Pi = Fi- m3V3U>2 + m2V2UJ3 

Newton : \ P2 = ^2 + "msVsUJi — niiViUJs , 
p^= F-i- m2V2UJi + miviuj2 

{vr^ = Ti + (m2 - m^)v2V-i + {h - h)^2i^2, 
■k^=T2 + (ms - mi)viV2, + {h - Ii)ojioj2, 
vTg = T3 + (mi - m2)viV2 + (/i - /2)a;ia;2 

showing brain's individual mass and inertia couplings. 



(8) 



^^In reality, is a 3 x 3 attitude matrix. However, for simplicity, we will stick to the (mostly) 
symmetrical translation-rotation vector form. 
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Equations ©-(IS 



energy of the brai: 



or, in tensor form 



Jl8 



can be derived from the translational + rotational kinetic 

(9) 



Ek = -v*Mv + -uj^Ilj, 



For this we use the Kirchhojf-Lagrangian equations: 
d 



dt 
d_ 

di 



d^Ek = d^Ek X cj + F, 
du>Ek = di^Ek -xu + dvEk x v + T, 



(10) 



where d^Ek = du^Ek = ^j-', in tensor form these equations read 



dt 
d_ 

di 



d,,E = elid,,E)u;'' + Fi, 

d.^E = El{d^,E)u^ + El,{d,,E)v^ + T,. 



Using (j9])- (|10p . brain's Unear and angular momentum co- vectors are defined as 

p = 5v-Efc, TT = dujEk, 

or, in tensor form 

Pi = d^rE, -Ki = d^rE, 
with their corresponding time derivatives, in vector form 



- =-d E 
dt dt ^ ' 



or, in tensor form 



d 



dt 



dyiE, 



vr,- =— TT,- = —o^iE, 



~dt ' dt 



^*In a fully-coupled Newton-Euler brain dynamics, instead of equation we would have brain's 
kinetic energy defined by the inner product: 



Eh 



\M 



SE(3) 
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or, in scalar form 



P= [Pi,P2,Pz] = [miVi,m2V2,m'i,V'i\, 77= [71-1,712,71-3] = [Iiioi, 12(^2, h'^z]- 

While brain's healthy Sii^(3)— dynamics within the cerebrospinal fluid is given 
by the coupled Newton-Euler micro-dynamics, the TBI is actually caused by the 
sharp and discontinuous change in this natural SE{'i) micro-dynamics, in the form 
of the S'i?(3)— jolt, causing brain's discontinuous deformations. 



6.2.4 Brain's traumatic dynamics: the 5£'(3)— jolt 

The ^-^(S)— jolt, the actual cause of the TBI (in the form of the brain's plastic 
deformations), is defined as a coupled Newton+Euler jolt; in (co)vector form the 
5^(3) -jolt readffl 



5^(3) - jolt : 



Newton jolt :F = p — pxcj — pxcj, 

Euler jolt : T = -ff — tt xcj — ttxcL; — pxv — pxv. 



where the linear and angular jolt co- vectors are 

F = Mv = [A, F2, hi t = Icj = [Ti, Ta, is], 

where 

are linear and angular jerk vectors. 

In tensor form, the 5-E'(3)— jolt readj^ 

Fi = Pi- elkPj^'' - ^kPi'^^^ 3,k = l,2, 3) 

fi = -el^i^ju^ -el^T^jUj^ -el^pjv^ -el^pjv\ 

in which the linear and angular jolt covectors are defined as 

F = Fi = m^ = Mi,v^ = [Fi,F2,F3], 

T = f, = Iu,=l,jO^ = [ti,f2,t3], 

where v = i;*, and u = uf" are linear and angular jerk vectors. 



^^Note that the derivative of the cross-product of two vectors follows the standard calculus 
product-rule: J^(uxv) = uxv + uxv. 

^"in this paragraph the overdots actually denote the absolute Bianchi (covariant) time-derivative 
((6]), so that the jolts retain the proper covector character, which would be lost if ordinary time 
derivatives are used. However, for the sake of simplicity and wider readability, we stick to the same 
overdot notation. 
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In scalar form, the ^^^(S)— jolt expands as 



Newton jolt : < 



Fi = pi - m20J3V2 + ms {0J2V3 + V3i02) - m2V20J3, 

F2 = P2 + miujsvi — msLOiVs — m^vsoji + miVidj^, 
^ -F3 = P3 — miOJ2Vi + m2UJiV2 — V2OJ1 — miviuj2, 



Euler jolt : < 



Ti = 7f 1 - (m2 - ma) {V3V2 + V2V3) - (h - h) (w3cj2 + ^2^3) , 
^2 = 7f2 + (mi - ms) (?;3-ui + viVs) + (/i - I3) {uj^uji + uiUs) , 
^3 = 7^3 — (mi — m2) {v2Vi + V1V2) — {h — h) (^2^1 + a;ia;2) • 



Wc remark here that the linear and angular momenta (p, tt), forces (F.T) and 
jolts (F,T) are co-vectors (row vectors), while the linear and angular velocities 
(v,a;), accelerations (v, w) and jerks (v, a>) are vectors (column vectors). This 
bio-physically means that the 'jerk' vector should not be confused with the 'jolt' co- 
vector. For example, the 'jerk' means shaking the head's own mass-inertia matrices 
(mainly in the atlanto-occipital and atlanto-axial joints), while the 'jolt'means ac- 
tually hitting the head with some external mass-inertia matrices included in the 
'hitting' SE(3)-jolt, or hitting some external static/massive body with the head 
(e.g., the ground - gravitational effect, or the wall - inertial effect). Consequently, 
the mass-less 'jerk' vector represents a (translational+rotational) non-collision effect 
that can cause only weaker brain injuries, while the inertial 'jolt' co-vector represents 
a (translational-|-rotational) collision effect that can cause hard brain injuries. 



6.2.5 Brain's dislocations and disclinations caused by the SE{3)—io\t 

Recall from introduction that for mild TBI, the best injury predictor is considered 
to be the product of brain's strain and strain rate, which is the standard isotropic 
viscoelastic continuum concept. To improve this standard concept, in this subsec- 
tion, we consider human brain as a 3D anisotropic multipolar Cosserat viscoelastic 
continuum, exhibiting coupled-stress-strain elastic properties. This non-standard 
continuum model is suitable for analyzing plastic (irreversible) deformations and 
fracture mechanics in multi-layered materials with microstructure (in which slips 
and bending of layers introduces additional degrees of freedom, non-existent in the 
standard continuum models. 

The S'£^(3)— jolt (F, T) causes two types of brain's rapid discontinuous deforma- 
tions: 

1. The Newton jolt F can cause micro-translational dislocations, or discontinuities 
in the Cosserat translations; 

2. The Euler jolt T can cause micro-rotational disclinations, or discontinuities in 
the Cosserat rotations. 
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To precisely define brain's dislocations and disclinations, caused by the ^-^(S)— jolt 
(F,T), we first define the coordinate co-frame, i.e., the set of basis 1-forms {dx*}, 
given in local coordinates X — (x ^ X ^ X ) = {x,y,z), attached to brain's center- 
of-mass. Then, in the coordinate co-frame {dx^} we introduce the following set of 
brain's plastic-deformation-related S" £'(3)— based differential p— formj^: 

the dislocation current 1-form, J = Jidx^; 
the dislocation density 2-form, a = ^aij dx"^ A dx^; 
the disclination current 2-form, S = ^Sij dx* A dx^; and 
the disclination density 3-form, Q = ^Qijk dx^ A dx^ A dx^ , 
where A denotes the exterior wedge-product. These four S'£'(3)— based differen- 
tial forms satisfy the following set of continuity equations: 

a = -dJ - S, (11) 

Q = -dS, (12) 

da = Q, (13) 

dQ = 0, (14) 

where d denotes the exterior derivative. 

In components, the simplest, fourth equation (jl4p . representing the Bianchi iden- 

Differential p— forms are totally skew-symmetric covariant tensors, defined using the exterior 
wedge-product and exterior derivative. The proper definition of exterior derivative d for a p— form 
/? on a smooth manifold M, includes the Pomcare lemma: d{dj3) = 0, and validates the general 
Stokes formula 

I I dp, 

J BM JM 

where M is a p— dimensional manifold with a boundary and dM is its (p— 1)— dimensional boundary, 
while the integrals have appropriate dimensions. 

A p— form jS is called closed if its exterior derivative is equal to zero, 

d/3 = 0. 

From this condition one can see that the closed form (the kernel of the exterior derivative operator 
d) is conserved quantity. Therefore, closed p— forms possess certain invariant properties, physically 
corresponding to the conservation laws. 

A p— form /3 that is an exterior derivative of some (p — 1)— form ct, 

P = da, 

is called exact (the image of the exterior derivative operator d). By Poincare lemma, exact forms 
prove to be closed automatically, 

dp = d{da) = 0. 

This lemma is the foundation of the de Rham cohomology theory 
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tity, can be rewritten as 

dQ = diQ[ijk] dJ- A dx" A dx^ A dx^ = 0, 

where di = while 0[ij...] denotes the skew-symmetric part of ^ij...- 

Similarly, the third equation p3|) in components reads 

-^Qijkdx'' Adx^ Adx^ = dka^ij] dx'^ A dx"^ A dx^ , or 
Qijk = — 6c?fca[ij]- 

The second equation (jl2|) in components reads 

^(jjjA: c^a;* A cix-' A dx'' = -dkS^ij-^ dx^ A dx* A dx^ , or 
Qijk = 65fe5[jj]. 

Finally, the first equation (jlip in components reads 

1 . . 1 . . 

-djj dx* A dx^ = {djJi — i^Sij) dx^ A dx^ , or 

In words, we have: 

• The 2-form equation (jlip defines the time derivative a =\aij dx* A dx^ of the 
dislocation density ol as the (negative) sum of the disclination current S and 
the curl of the dislocation current J. 

• The 3-form equation (fT2]) states that the time derivative Q =-^Qijk dx'' Adx^ A 
dx^ of the disclination density Q is the (negative) divergence of the disclination 
current S. 

• The 3-form equation (fT3]l defines the disclination density Q as the divergence 
of the dislocation density ex, that is, Q is the exact 3-form. 

• The Bianchi identity (jl4p follows from equation (|13p by Poincare lemma and 
states that the disclination density Q is conserved quantity, that is, Q is the 
closed 3-form. Also, every 4-form in 3D space is zero. 

From these equations, we can derive two important conclusions: 
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1. Being the derivatives of the dislocations, brain's discHnations are higher-order 
tensors, and thus more complex quantities, which means that they present a 
higher risk for the severe TBI than dislocations — a fact which is supported 
by the literature (see review of existing TBI-models given in Introduction of 

m)- 

2. Brain's dislocations and disclinations are mutually coupled by the underlaying 
^-^(S)— group, which means that we cannot separately analyze translational 
and rotational TBIs — a fact which is not supported by the literature. 

For more medical details and references, see [25] . 
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